VLASOV-MAXWELL-BOLTZMANN DIFFUSIVE LIMIT 



JUHI JANG 



Abstract. We study the diffusive expansion for solutions around Maxwellian 
equilibrium and in a periodic box to the Vlasov-Maxwell-Boltzmann system, 
the most fundamental model for an ensemble of charged particles. Such an 
expansion yields a set of dissipative new macroscopic PDE's, the incompress- 
ible Vlasov-Navier-Stokes-Fourier system and its higher order corrections for 
describing a charged fluid, where the self-consistent electromagnetic field is 
present. The uniform estimate on the remainders is established via a unified 
nonlinear energy method and it guarantees the global in time validity of such 
an expansion up to any order. 



1. Introduction and Formulation 

The dynamics of charged dilute particles can be described by the celebrated 
Vlasov-Maxwell-Boltzmann system: 

d t F+ + v ■ V X F+ + (E + v x B) ■ V V F+ = Q{F+, F+) + Q{F+, F_), 

' ' ' ' ' d t F- + v ■ V X F- - (E + v X B) ■ V„F_ = Q(F_, F+) + Q(F_,F_), 

with initial data F±(0, x,v) — Fo ±(x,v). For notational simplicity we have set all 
physical constants to be unity, see [Hj for more background. Here F± (t, x, v) > are 
the spatially periodic number density functions for the ions (+) and electrons (-) 
respectively, at time t > 0, position x = (xi,X2,x^) e [ — it,it] 3 = T 3 and velocity 
v = (vi,V2,V3) € M 3 . The collision between particles is given by the standard 
Boltzmann collision operator Q{G\ 1 G2) with hard-sphere interaction: 

(1.2) Q(G 1 ,G 2 )= f \(u-v)-io\{G 1 {v')G 2 (u')-G 1 (v)G 2 {u)}dudu, 

JR 3 x S 2 

where v' = v — [(v — u) ■ lo]lo and v! = u + [(v — u) ■ u]w. 

The self-consistent, spatially periodic electromagnetic field [E(t, x), B(t, x)} in 
(1.1) is coupled with F±(t,x,v) through the Maxwell system: 



(1.3) 



d t E~V x B = - / v(F + - FJ)dv, V • B = 0, 
Jr 3 

d t B + S7 x E = 0, V-E= [ (F+-F-)dv, 

Jr 3 



with initial data ^(0, a;) = E (x), B(0, x) = B (x). 

It turns out that it is convenient to consider the sum and difference of F + and 
F_ as proposed in 0|. Defining 

(1.4) F = F + + F- and G = F + — F-, 
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(1.1) and (1.3) can be rewritten as following: 

d t F + v ■ V X F + (E + v x B) ■ V V G = Q(F, F), 
d t G + v- V X G + {E + v x B) ■ V V F = Q(G, F), 

(1-5) d t E-X7xB = -[ vGdv, V • B = 0, 



d t B + V x E = 0, V-E= / Gdv, 



with initial data F(0,x,v) = Fq(x,v), G(0,x,v) = Gq(x,v), E(0,x) = Eq(x) and 
B(0,x) = B Q (x). 

Now we introduce the diffusive scaling to (|1.5|l : for any e > 0, 

e5 t F £ + u • V x .F e + (E £ +vx B £ ) ■ V V G £ = -Q(F £ .F £ ), 

e 

ed t G e + v ■ V X G £ + {E e + v x B e ) ■ V V F £ = -Q(G £ ,F £ ), 

(1.6) r 

ed t E £ - V x B £ = - \ vG e dv, V • B e = 0, 

ea tJ B e + V x E £ = 0, V • £ e = / G e efa. 



For notational simplicity, we normalize the global Maxwellian as 
< L7 > - 12^'-"'"' 2 - 

We consider the following formal expansion in e around the equilibrium state 
[F, G, E, B] = 0, 0, 0]: for any n > 1, 

F £ (t, x,v) = [i + v7*{e/i(t, x, v) + e 2 f 2 (t, x,v) + ...+ e n f £ (t, x, «)}, 

G e (t, u) = ^{egi (t, x, u) + e 2 g 2 (t, x, v) + ... + e n g £ l (t, x, v)}, 
(1.8) „ 

E £ (t, x) = {eEi{t, x) + e 2 E 2 (t, x) + ... + e n E £ n {t, x)}, 

B £ {t,x) = {eB l {t 1 x)+e 2 B 2 {t,x) + ... + e n B £ n (t,x)}. 

To determine the coefficients /i(t, x, v), / n _i(t, x, u); <7i(t, x, u), g n ~i(t, x, v); 
Exit, x), E n _x(t, x); Bi(t,x), B n _i(t,x), we plug the formal diffusive expan- 
sion l|1.8[) into the rescaled equations (|1.6fl : 

(1.9) 

(eft + t ;-V :c ){£/ 1 + ...+£"/ r t} 

+ ^{e(Fx + v x B{) + ... + e n (Ef l +vx B £ n )} ■ V v [^J2{e gi + ... + e n g £ n }] 

' Q(n + ^{e.h + ... + s n f £ }, n + ^{eh + ... + e n f*}), 



(sd t + vV x ){egi + ... + s n g s n } 

+ ^ F {e(E 1 +vxB 1 ) + ... + e n (E £ n +vx B £ n )} ■ V v [fx + tfi{efi + ... + e n f £ }] 
V f 1 

1 Q(^{e.9i + ... + e n g £ n }, /x + VT^h + - + e n /»}), 
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sd t {eE 1 + ... + e n E £ n } - V x {eB 1 + ... + e n B e n } = - [ v^J2{e gi + ...e n g e n }dv, 

J~R 3 

ed t {eB 1 + ... + e n B e n } + V x {eE l + ... + e n E e n ) = 0, 

V • {eE 1 + ... + S n E s J = / ^{e 9l + ...s n g e n }dv, 

V ■{eB 1 + ... + E n B s n } = 0. 

To expand the right hand side Q in the above, we define L the well-known linearized 
collision operator and C another linearized operater as 

(1.10) Lf = --^={Q(n,^lf) + Q(^Jlf,n)}, 

(1.11) £ 5 = --^Q(Vm 5 ,m), 

and the nonlinear collision operator T as (non-symmetric) 

(1-12) r(f,g) = -j=Q(y/jif,y/jig). 

Note that Lf and Cg can be written as following in terms of T: 
(1-13) L/ = -{r(VM,/)+r(/,VM)}, Cg = -T(g,^Jl). 

Now we equate the coefficients on both sides of the equation (11.911 in front of 
different powers of the parameter e. Let 

/_i = /o=0, <7_i=<?o = 0, E = 0, B = 

to obtain 

d t fm + v ■ V x f m+ i + -j= ^2 ( E i + vx B 'i) ■ V v {«/jigj) 

= -Lf m+2 + r (/*>/i)' 
(1-14) x 

9t5m + U ■ Vagm+l + — (Ej + V X Si) ■ V v (y/fifj) - E m+ i ■ V^/JI, 

V " i+j= m +l 
i,i>l 

= -£.9m+2 + X! r (9hfj), 



i+j— m+2 

for — 1 < m < n — 3 as well as 



(1.15) 



9 t -B m - V x B m+1 = - J vg m+ly ^dv, V ■ B m+ i = 0, 
<9 t 5 m + V x E m+1 = 0, V • £ m+ i = J g m+1 ^/JIdv, 



for < to < n — 2. Moreover, we can collect terms left in <|1.9ll with powers e n 1 or 
higher and divide by e™ _1 to get the equations for the remainders ^^g s n ,E^,B e n . 



4 



JUHI JANG 



Note that all the e m+1 -th order terms vanish for m < n~ 3 because of (|1.14|) . First, 
we write equations for and g e n : 

(1.16) 

e 2 d t f n + ev ■ V x f* + Lfn = 



{-dtfn-2 - v ■ VJn-i - j= (Ei + vxB i )-V v ( y /Jlg j )+ ]T Tfafj)} 

i-\-j=n—l 1+3=71 

> 1 

+ e{-d t f n -i -yz J2(Ei + vx Bi) ■ V^(VMft) + E &)} 

i,j>l 

n-1 

+ e T(/«, /») + ]T e'{r(/«, /<) + r(/i, /*)} + X! ^ +J ~" r (/» /i) 

i=l i+j>n+2 

£ n+l 

-—(E s n +vxBl)-V v (,/jlg s n ) 

^ n— 1 

-77 E ^ n+1 (£,+«B,)-v„(^); 

v " i+i>n+l 



£ <9 t c£ + • V x g„ - eE e n ■ v^fji + Cg e n = 

{-9 t g n -2 - V ■ VxSn-l -= ^ (Ei+V X Si) -Vvii/jlfj) + E n -i -Vy/Jl 

i,i>l 

+ £ 1%, /,•)} + £{-9t3n-l " y= E + « X ^) • v^Vm/j-) 

2+_7=fi V ' i-\-j—n 

i,j~> 1 i ?< 7>l 

+ 2 r( 9i) / J )} + e n r(^ s /=) + 53 e < {r( ff =,/ i ) + r( 9i) /«)} 

i+j— n+1 z— 1 

t,i>i 

+ £ ^-Tfe./j) - — + « x B^) ■ 

i+j">n+2 V r 1 

1 n — 1 

- -f E £l+1 {(^ + « x Bi) ■ V„(vft£) + + v x fl£) • V„(V7^)} 
* »+3>n+l 

Similarly, by (|1.15(1 . we get the remainders for E^, B e n : 

ed t El - V x B s n = - d t E n - X - [ v&Vfidv, V • B e n = 0, 
(1.17) Jr3 

sd t B e n + VxE E n =-d t B n _ 1 , V-E e n = / g s n y/JIdv. 

Jm 3 
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The fluid equations can be obtained through the conditions l|1.14J) and Ijl.l5|l . 
We first recall that the operator L > 0, and for any fixed (t, x), the null space of L is 
generated by [y/j2, v^/JL, |«| %//!]. For any function f(t, x, v) we thus can decompose 

/ = P 1 / + {I-P 1 }/ 

where Pi/ (hydrodynamic part) is the L? v projection on the null space for L for 
given (t, x). We can further denote 

(1.18) Pi/ = { Pf (t,x)+v u f (t,x) + (^- - l)0 f (t,x)}^jl. 
Here we define the hydrodynamic field of / to be 

[p f (t,x),Uf(t, x),8 f (t,x)} 

which represents the density, velocity and temperature fluctuations physically. For 
the velocity field u / (t, x) , we further define its divergent-free part as 

Pquj = the divergent-free projection of Uf 

so that 

V • {P u f } = 0. 

Similarly, one can show that C > and for any fixed (t,x), the null space of C 
is one dimensional vector space generated by Likewise, any g(t,x,v) can be 

decomposed into 

.g = P 2 .g + {I-P 2 } 5 

where P23 (hydrodynamic part) is the L\ projection on the null space for C for 
given (t, x). We can further denote 

(1.19) P 2 g = o- g {t 7 x)^JJ. 

Here o g {t, x), the hydrodynamic field of g, can be interpreted as the concentration 
difference. For more details about C and P 2 , we refer [3]. Before going on, we state 
the coercivity of L and C which will be often used in the subsequent sections: there 
exists a 8 > such that 

(1.20) (Lf,f) ><5|(I-Pi)/|2, {Cg,g)>S\(I-P 2 )g\l. 

See Lemma 1 in [S] for its proof. Note that the operator L defined in 8 is equivalent 
to [L, C] in our case. 

Now define [p m , u m , 6 m , a m ] to be the corresponding hydrodynamic field of the 
m-th coefficients f m and g m . As for the first coefficients fi(t,x,v) and gi(t,x,v), 
from (|1.14|) 

(1.21) {I - Pi}/i = and {I - P 2 }. 9l = 
which immediately yield that 

(1.22) B x = and E x = V0i 

up to constant and for some function <f>i(t,x) satisfying A</>i = <ti; in particular, 
B\ may be assumed to be zero physically in a sense that nonzero constants B\ do 
not cause the hydrodynamic equations H1.25|) - (|1.28|) to change. It will be shown in 
Lemma \'6. II that its velocity fluctuation Ui(t, x) is incompressible: 

(1.23) V • ui = or ui = P ui, 
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and its density and temperature fluctuations pi{t, x) and 0i(i, x) satisfy the Boussi- 
nesq relation: 

(1.24) Pl +e 1 = o. 

Moreover, [ui,0i,oi] satisfies the nonlinear incompressible Vlasov-Navier-Stokes- 
Fourier equations: 

(1.25) d t U\ + Ux ■ Viti + Vpi = rjAui + 0iV<^>i, 

(1.26) d t ui + i*i • Voi = aAcri — aa\, 

(1.27) Acf) 1 =a 1 , 

(1.28) df0i+ui-V0i = «A0i, 

where pi(t,x) is the pressure and 77, k, a > are physical constants. 

As for the coefficients f m {t, x, v), g m (t, x, v) for m > 2, by l|1.14|l . the microscopic 
part of f m and g m is determined by: 

{I - P x }/ m - L-M-Qt/m-a - « • V x / m _! + ^ r(/i, /,) 

(1.29) -- L £ (£ !+u xB 1 ).V„(^ ft )K 
{I - P2}.9™ = £~ 1 {-d t g m -2 - v ■ V x g m -i + ^ T ^9iJj) 

i+j=m 
i,J>l 

(1.30) ^ (£ 1 +«B 1 ).V„(V)]/ J )+£ m -i-^}. 

i,i>l 

On the other hand, for the hydrodynamic field of / m and g m : 

M 2 3 

Pl/m = {pm(t,x) + V ■ U m (t,x) + {— -}9 m (t,x)}^/p, 

~P?g m = a- m (t,x)y/JI, 

we can deduce an m-th order incompressibility condition 
(1-31) V • {I- Po}u m = -dtp m -l, 

an m-th order Boussinesq relation 
(1.32) 

Prn + 0m = A _1 V ■ {-U t ■ V(P U ro _l) - P()U m -l ' VUl + E\O m -X + E m _ x O x 
+ Kn-l} + (^^, -^({I - Pl}« ' VxPi/m-l)) - - • Ul, 

and an m-th order linear Vlasov-Navier-Stokes-Fourier system for [Pou mi m , a m ]: 

(1.33) (d t + ui ■ V - 7?A)P Wm + foUm ■ Vtti + V Pm - [E\a m + E m a x ) = R u m , 

(1.34) (ft + ui • V - aA + a)cr m + P u m • Vox = iC, 

(1.35) (ft +ui • V - nA)9 m + P u m • V0 X = i? m , 
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(1.36) V X E m = -dtB m -i, V-E m =a m , 

(1.37) VxB m = [ {I-P 2 } g m v^fldv + d t E m -i, V • B m = 0, 

JR 3 

with compatibility conditions coming from conservation laws 

(1.38) / E m dx — -a / E m dx + £ m -i, [ B m dx = 0. 

Here i?^, i?^, i?^ and defined precisely in l|3.2[) . I|3.4|l and l|3.20|l . essen- 

tially depend only on /,,-, g h E } , Bj for j < m- 1, since {I - Pi}/ m , {I - P2}ffm, 
{I — Po}u m , as well as p m + 9 m have been determined. 

In order to state our results precisely in the next section, we introduce the 
following norms and notations. We use (• , ■) to denote the standard L 2 inner 
product in R 3 , while we use (• , •) to denote L 2 inner product either in T 3 x M 3 
or in T 3 with corresponding L 2 norm || ■ ||. We use the standard notation H s to 
denote the Sobolev space W s ' 2 . For the Boltzmann collision operator 11.2f> . define 
the collision frequency to be 

(1.39) v(v) = / \v - v'\n(v')dv', 

JR 3 

which behaves like \v\ as \v\ — ► 00. It is natural to define the following weighted L 2 
norm to characterize the dissipation rate. 

2 — / „2L,\,,/„,\J», ll„l|2 — / „2 



\g\Z= / g 2 {v)v(v)dv, \\g\\i= / g 2 {x,v)v{v)dvdx. 

JR 3 JT 3 xK 3 

Observe that for hard sphere interaction, 

(1-40) \\(l + \v\)ig\\<C\\g\\„. 

In order to be consistent with the hydrodynamic equations, we define 

(141) d p = d 11 d 12 d 13 d Pl d p2 d Pa 

where 7 = [71,72,73] is related to the space derivatives, while [3 — [/3i, /?2, /%] is 
related to the velocity derivatives. 

We now define instant energy functionals and the dissipation rate. 

Definition 1 (Instant Energy) For N > 8, for some constant C > 0, an 
instant energy functional £^(f,g,E,B)(t) = £jsi(t) satisfies: 

(1.42) if w (t) < Y, H^/> ^9]\?{t) + J2 \\[SyE, d,B}\\ 2 (t) < CS N (t). 

\P\+\l\<N h\<N 

Definition 2 (Dissipation Rate) For N > 8, the dissipation rate T>jsr(f,g)(t) 
is defined as 

v N (f,g)(t)= Y III^Pi/: 5 7 P 25 ]|| 2 W 

| 7 |<JV 

(1.43) j 

+p» E ll[^{I-Pi}/, d p {I-P 2 }g]\\l(t). 

\0\ + h\<N 

We remark that both the instant energy and the dissipation rate are carefully de- 
signed to capture the structure of the rescaled Vlasov-Maxwell-Boltzmann equation 
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(|1.6(l . First of all, the electromagnetic field [E, B] is included only in the instant 
energy, which prevents the exponential decay on £n unlike the pure Boltzmann 
case for hard potentials. See |U] and [T2|. Notice that there is no ^ factor in front 
of the hydrodynamic part [Pi/, P25] m the dissipation rate 2?jv (/, 3), since only 
the microscopic part [{I — Pi}/, {I — P2j-g] should vanish as e — > 0. For notational 
simplicity, the Einstein's summation convention is used for Greek letter up to order 
N > 8 throughout the paper, unless otherwise specified. We denote V = W x and 
use C to denote a constant independent of e. We also use U(-) to denote a general 
positive polynomial with U (0) = 0. 



2. Main Results 



The first result is to determine the coefficients /1, /2, ■■■fm] <?i, <?2, ■■■9m] E\, E2, 
■ ■■E m ; B\,B2, ■■■B m in a diffusive approximation (|1.8|) . 

Theorem 2.1. Let m divergent-free vector-valued functions \u\(x), u^x), u^ n (x)], 
2m scalar functions [9i(x), 02(x)...0^ n (x); a® (x), o^ix), ...cr„(a;)] be given such that 



and 



lKlliP + l|0illi^ + lkillff3 < a/ 

/ CT-°(:z;)dx = 0, / u° r (x)dx = - V / E°(x) x B°(x) dx, 
Jt 3 Jt 3 ij~t r h 3 

/or 1 < r < m. Here E® , B^ 1 are defined by Ei(0, x), Bj(0, x) which have been induc- 
tively determined at the precedents since i,j < r, starting with the average condi- 
tions J T 3 Bidx = and J T3 E^dx — 0. T/ien /or sufficiently small M and given m 
real vectors e±(= 0), e2, e m , t/iere exist unique functions fi(t, x, v), fi{ti x, v), 
f m (t,x,v); gi(t,x,v),g 2 (t,x,v),...,g m (t,x,v); Ex(t,x),E 2 (t,x), ...,E m (t,x) and 
Bi(t,x),B 2 (t,x), ...,B m (t,x) with 

/ ay(t, x)dx = 0, / Pou r {t 1 x)dx — — 2, / Ei(t,x) x Bj(t,x) dx, 
Jt 3 Jt 3 

, . i . 7 > 1 

(2.1) 

/" 3 1 x /" 

/ -9 r {t,x)dx = -- / Ei(t,x) ■ Ej(t,x) + Bi(t,x) ■ Bj(t,x) dx, 

Jt 3 2 2 .■ t~„ Jt 3 



i+j=r 



suc/i iftai initially Pou r (0,x) — u®(x), 9 r (0,x) — 0^(x), a r (0,x) = <?r( x ) and 

f T 3 E r (0, x)dx = e r , and fi(t, x, v), gi(t, x, v), E\(t, x), Bi(t, x) satisfy 

and f r (t,x,v),g r (t,x,v),E r (t,x),B r (t,x) satisfy \1.29^ - \1.3t$i for 2 < r < m. 
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Moreover, for 1 < r < m, for any (3, and for all s > 3, there exists a polyno- 
mial U r ,f3, s with U r ,f3,s{0) = such that 



J2{\\[9?fr, d?g r Mt) + \\[d T E r , d T B r }\\(t)} 



(2.2) M ^ S 

<e A 'C/ r ,/3,s( ^ {||Wj||j?2»+4(r-i) + ||0°||jJ2»+4(r-J) + ||<7-°||jj2«+4(r-i)}), 
l<j<r 

where space-time derivatives 

and A can be chosen as | min{?7, k, a} for sufficiently small M. 

We now turn to the most important question about the remainder estimates for 
f £ , g £ n , E 6 , and B^. We first study the classical case for the first order remainders 

f e = fl, g £ = gf, E £ = Ef, B £ = Bf 

which satisfy the nonlinear Boltzmann type equations: 

(2.3) d t f + -v ■ V x f + \hf = -T(f, f) - -^-(E £ + vxB £ )- V v (^J2g s ), 

£ £ £ \f i-L 



i ,_ , 1 _ i„, „ _ l 



8 t g £ + -v ■ (V x g £ - ^/JlE £ ) + -Cg £ = -T(g £ , f £ ) -—(E £ +vxB £ )- V^f), 



ed t E £ - VxS £ = - [ g £ Vy/JI dv, V • B e = 0, 
(2-4) J 

ed t B e |VxE E = 0, V • E £ = / g £ ^/JJ dv. 

Theorem 2.2. Let N > 8. Let f £ (0,x,v) = f$(x,v), g £ {0,x,v) = g £ (x,v) 
E £ (0,x) = E £ (x), B £ (0,x) — B e (x) satisfy the mass, momentum and energy con 
servation laws: 

(/ £ ,Vm) = o, (sg,v^) = o, 

^ 2 5 ) (/o) v \^) + 6 J E o x B o dx = 0, 

(/oMH 2 Vm)+^ / \E £ \ 2 + \B £ \ 2 dx = 0. 

Then there exists an instant energy functional £n{S £ ,g £ \E £ ,B £ )(t) such that 
if £N(f £ ,g £ ,E £ ,B £ )(0) is sufficiently small, then 

(2.6) J t £N{r ' E ^ B£){t) + Pw(/£ ' ff£)(<) - °' 
In particular, 

(2.7) sup E N {f £ ,g £ ,E £ ,B £ )(t)<£ N {f £ ,g £ ,E £ ,B £ ){Q). 

0<*<oo 

Moreover, for k > 1 t/iere exists Cjv.fc > such that 

(2.8) £ N (f,g E , E £ , B £ )(t) < C N , k £ N+k {f, g £ , E £ , B £ )(0) |l + ^ 
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We remark that our initial data /g , g$ , Eq , Bq are general and can contain ini- 
tial layer for the Vlasov-Navier-Stokes-Fourier limit. We can easily include (one) 
temporal derivative dt in our definition of instant energy and dissipation rate, and 
obtain the same uniform bound for such a new norm. With such a modification, 
the boundedness of 9t/g, dtgo, dtEg, dtB^ automatically removes the formation 
of any initial layer. 

For higher order remainders g e n , E^ 7 B^ with n > 2, we have 

Theorem 2.3. Let N > 8. Given /i,/ 2 ,.../„; 9i,92,-9n, E t E 2 ,...E n ; B x B 2 ....B n 
constructed in Theorem \2.1\ and let 

Ill[ u n>0n! cr n]|lliv(*) 

(2.9) _ {\\u°\\ H 2N + 10+Hn-j) + \\Oi\\ H 2N + 10+Hn-j) + | |of | | #2JV + 10+4( Il -j) } . 

l<j<n 

And let 

F e (0, x, v) = fi + VM{e/i(0, x, v) + ... + e—Vn-iCC x, «) + ^(0, z, «)}, 
G £ (0,a;,u) = ^{e<7i(0, x, i;) + ... + e™" V-i(0, a, u) +^(0,1,8)}, 
£ £ (0, x) ee eSi(0, a:) + ... + e n - 1 £^_ 1 (0, a;) + e n ££(0, x), 
B E (0, x) ee eB^O, x) + ... + e"- 1 B n _ 1 (0, a) + e"B* (0, i) 
be given initial data satisfying the following conservation laws: 

{F s (0,x,v) - fi(v)}dvdx = 0, / / G e (0, x, v)dvdx = 0, 
1 Jt 3 Jr 3 



(2.11) / / vF e (0,x,v)dvdx + E £ (0, x) x B e (0, x)dx = 0, 
Jt 3 Jr 3 jt 3 

|«| 2 {F E (0,a;,u)-/x(«)}cfod2; + / |£ e (0, :r)| 2 + |.B £ (0, x)| 2 dir = 0. 

Jt 3 

Then there exist an instant energy functional £n and a positive polynomial U with 
U(0) = such that if both e and 

<^v(/n - fn, 9n ~ 9n, E„ - E n , B^ - B n )(0) 

are sufficiently small, then 

SUp £ N (fn ~ fn, 9„ - 9n, E„ - E n , B^ - B n ){t) 
0<t<oo 

(2.12) < {e u m ul MmN)EN{rn _ fn gl _ g ^ K _ K _ Bn){Q) 

+ e 2 U(\\\[u n ,e n ,a° n }\\\ N )}. 
Moreover, for k > 1 there exists Cjv,fc > such that 

&N{fn ~ fn, 9n ~ 9n, E^ - E n , B £ n - B n )(t) 

(2 13) < C N>k {e u ^<^^MS^ +k (f^ - fn,9 e n -9n,E e n - E n , B e n - B n )(0) 
+ £ 2 {/(|||[u°,0°, ( T°]||| Ar )}|l + ^ 

Note that the conservation laws (|2.11|l in Theorem 12.31 imolv the conservation 
laws for g e n , E^, B^ due to (|2.1(l . See l|7.5(l for the precise description. In 
addition, we remark that the positivity for the initial data (F £ ± G £ )(0,x,v) > 
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which is equivalent to F±(0, x, v) > through the relation ((1.4(1 can be verified in 
the same way as discussed in the Appendix of [H]- 

There are several macroscopic fluid models for classifying the dynamics of a 
charged fluid, but none has been derived from the Boltzmann theory mathemat- 
ically. This is because the construction of the global solution to the important 
Vlasov-Maxwell-Boltzmann system had been open for a long time until only a few 
years ago, in jS], a unique global in time classical solution near a global Maxwellian 
for such a master system was constructed. However, it still remains a major open 
problem: to construct global renormalized solutions to the same system. 

Singular limit problems emanating from the Boltzmann equations have been 
studied by many people for decades [H El El El El El El EH EI 1121 [HI ESI El- 
a great overview of the issue is given in |§1 1181 IT%| . In particular, in the recent 
work [5], higher order approximations with the unified energy method have now 
been shown to give rise to a rigorous passage from the Boltzmann equation to the 
Navier-Stokes-Fourier systems beyond the Navier-Stokes approximation. 

In this article, we rigorously establish the global in time validity of the diffusive 
expansion l|1.8|l to the rescaled Vlasov-Maxwell-Boltzmann equations 6(1 for any 
order. This not only gives the uniform estimates for the higher order corrections, 
but also leads to the mathematical derivation of new dissipative hydrodynamic 
equations, which we call Vlasov-Navier-Stokes-Fourier System. An interesting fea- 
ture of our result is that for the chosen incompressible regime, the nontrivial mag- 
netic effect in the limit first occurs in the second order expansion (n — 2), while the 
electric effect is important at any order. We believe our result opens a new line of 
research for those macroscopic approximations both physically and mathematically. 

The method of this paper is based on the improvement of the recently developed 
nonlinear energy method in [HI El- We use the reformulation 1)1. 5(1 . by introduc- 
ing new unknowns 11.4(1 . to simplify both the character of the Vlasov-Maxwell- 
Boltzmann system and the analysis presented in this paper endowed with the clear, 
concrete linearized collision operators L and £. Macroscopic equations and local 
conservation laws are proven to be key tools to build the crucial positivity of L and 
C The most important analytical difficulty lies in deriving the uniform estimates 
on f £ and g e . Due to the singular behavior of the time derivative in our prob- 
lem, the positivity estimate for purely spatial derivatives is invoked: see 1(5. 4|l in 
Lemma f5.ll Notice that an additional term ^G(i) is needed. Such a differential 
form can still yield decay estimates with the notion of equivalent instant energies. 
The local conservation laws are used to estimate the more singular temporal de- 
rivative for the hydrodynamic field directly in terms of purely spatial derivatives 
of the microscopic part. We also use the trick of the integration by parts in the 
time variable and, in turn, by using such a uniform estimate, avoid encountering 
[9t{I — Pi}/ e , dt{l — P2}.9 £ ]- This complication is another reason for introducing 
the notion of equivalent instant energies. 

As pointed out in jS], the Vlasov-Maxwell-Boltzmann system is tricky to handle 
due to the hyperbolic nature of Maxwell equation and indeed, it is the most intrigu- 
ing and critical part of this article to control electromagnetic fields. For that, we 
first utilize the macroscopic equation 1(5.13(1 to estimate electric fields and then the 
Maxwell equation itself for magnetic fields. Noting that the macroscopic equation 
((5.13(1 is simpler than the one considered in 8 owing to the reformulation 1(1.5(1 . 
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we remark that our method together with the positivity estimate Ij5.4|l provides 
another lucid and concise way of proving the global in time classical solution for 
the Vlasov-Maxwell-Boltzmann equations, especially without using the temporal 
derivatives. 

On the other hand, it is delicate to establish the well-posedness of the new 
hydrodynamic equations because of their complexity, mainly stemming from elec- 
tromagnetic fields. It turns out that the compatibility conditions H1.38[l for aver- 
ages of electromagnetic fields are necessary and sufficient conditions in order for 
the Vlasov-Navier-Stokes-Fourier systems to have the unique solution. We should 
point out that those conditions are natural restrictions in that they can be derived 
from conservation laws. In turn, the exponential decay of hydrodynamic variables 
as well as the solvability of hydrodynamic equations lead to the solvability of kinetic 
equations and almost the same decay for approximate solutions. While the same 
exponential decay rate is obtained for the pure Boltzmann case, we are able to 
obtain only the polynomial decay rate l|2.8[l and l|2.13|l for solutions of the Vlasov- 
Maxwell-Boltzmann equations (compare with (2.6) and (2.13) in [2]) by applying 
the method proposed in [16| . In particular, for higher order remainders the more so- 
phisticated continuity argument is employed in order to get the desired polynomial 
decay rate l|2.13[l from l|7.15[> . 

The paper will proceed as follows: we will derive the high order linear Vlasov- 
Navier-Stokes-Fourier system in Section 3; we will prove Theorem l2.1l for coefficients 
fi(t,x,v),...,f m (t,x,v), gi(t,x,v),...,g m (t,x,v), Ei(t,x),...,E m (t,x), Bi(t,x),.., 
B m {t, x) in Section 4; Section 5 will be devoted to the positivity of L and C; in the 
last two sections, the first and higher order remainder estimates-Theorem 12.21 and 
Theorem 12 .31 - will be proven respectively. 



3. High Order Vlasov-Navier-Stokes-Fourier System: 
Formal Derivation 

In this section, we derive the microscopic equations (ll.29|) . I|1.30[l and hydrody- 
namic equations (|1.31() - (|1.38() . We shall use many results from 

Lemma 3.1. Assume that the expansion HI. 8]) satisfies H 1.6]) and such that 
for |t| + |/3| < N, 

n-l 

J2{\\[d?f m ,d^ 9m \\U+\\[9rE m ,d T B m ]\\} 

K 6 - 1 ) m=l 

+ 1 1 1 1, + 1 1 [drK, drB S n ] || < 00. 

Then there exist f n (t,x, v), g n (t,x,v ), f n+ i(t,x,v), g n+1 (t,x,v) and E n (t,x), 
B n (t,x) such that \1. 1J$ and 11.15]) are valid for all m < n — 1. Moreover, f\ 
and gi satisfy ]1.21]) . the incompressibility condition ]1.2'd\) . the Boussinesq rela- 
tion \1.2$ , an d the first order Vlasov-Navier-Stokes-Fourier system ]1.25]) - ]1.2ti]) : 
E\ and Bi satisfy \1.22\i . Form > 2, f m , g m , E m , B m satisfy the microscopic equa- 
tion and \1.H0\) . the m-th order incompressibility condition ]1.S1)) . the m-th 
order Boussinesq relation 11.32]) . and the m-th order Vlasov-Navier-Stokes-Fourier 
system f03)) - iT03)) with 
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R u m = {v- V x L~\d t {l - Pi}/ m _i - E T ^ fj)h v Vfl 

i+j=m+l, 

+ (V • V X £ _1 {{I - P l}(« ' V 4I - Pl}/m), «VM> 

- (v ■ VxL-Hrai, {i - Pi}/ m ) + r({i - h}/™, A)}, 

+ (v-V x L- 1 {-^= J! (£.+«B 1 ).V„(.9 3 ^)},"Vm) 

V i-\-j—m, 
i,j>l 

- (d t + ui • V - »yA){J - F }«m - {/ - Po}«m • Vui 

- (V ■ {/ - P a }u m ) Ul + |V(V • {J - P }u m ) 

+ Yl {EiV-E J -(d t E j - 1 -'VxB j )xB i }; 

iC = (« • V^-^ftU - P 2 }g m -i - E T (9i, fi)h VF) 

i+j=m+l, 

+ (v ■ V^-^jl - P 2 }(v • V X {I - P 2 } ffm ), VM> 

- (« • v x £ _1 {r(pi, {i - p 2 }.g m ) + r({i - p 2 }. 9m , 9l )}, Vm) 

+ <«■ V^" 1 ^ ^ (Bi+^xBi).^^)},^) 
V " iJ r j =m , 

- {V • (7 - P )Mm}<ri - (I - P )u m • Vai; 



j^^^.v^L-H^i-Pi}/™-!- E r(/ <5 /,•)}, 

+ (t, • V X L-\{1 - Px}( W • V X {I - Pi}/ m ), 

- (v ■ v^-Hr^i, {i - Pi}/ m ) + r({i - Pi}/™, A)}, J^^) 

i,3>l 

2 

+ (-d t - ui • V + KiA){p m + 6< m } 

- {V • (7 - Po)w m }#i - (J - Po)"™ • V01 
-\ E ^-(Pm-VxB^; 

i+j=m+l, 
»,J>1 
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(3.5) p m = ( Pm+1 + 6 m+1 ) - (Mj^, - PiXw • V a Pi/ m )}) 

5 



>U,so i/ie conservation laws \2. 1\) are valid for each 



Proof. Notice that l|1.14|l is clearly valid for m < n — 2 under the assumption (|3.1(l 
since it can be derived by letting e — ► after dividing l|1.9|l by e m+1 . By the same 
token, 11.15fl holds for m < n — 1. We now show the existence of the coefficients 
/„, g n , E n , B n and f n +i, 9n+i so that the equations (|1.14|) for m = n — 2, n — 1 and 
H1.15fl for m = n— 1 are satisfied. Firstly, by (|3.1|) . up to a subsequence, there exist 
/n, 9n, E n ,B n such that 

[/n.Pn] [fn,9n] weakly in || • \\» and BfJ [E n ,B n ] weakly in || • ||. 

The assumption H3.1|) guarantees the existence of their derivatives. By subtracting 
off Lf n , Cg n on both sides of (|1 . lf>|> we can isolate the zeroth order terms in 
the remainder equations and deduce (| 1 . 1 4|> for m = n — 2 at least in the sense of 
distributions. Moreover, letting e — > of l|1.17f) . we also deduce l|1.15|l for m = n— 1 
at least in the sense of distributions. 

Next, in order to find f n +i and g n +i, we use (|1.14(l for m — n — 2 to simplify 
(US) as 



I 



= -Edtfi ~ V ■ V x f n 



i-\-j=Ti i-\-j—n-\-l 

2,J>1 

(3.6) 

n-l 

+ e - 1 r(^, /*) + ^HWn, fd + m, /»)} + E e i+3 '- n_1 r(/<, /;) 

i=l i+j>n+2 

-— (^ + «X^).V (VM5n) 
^ n— 1 

- -F E ^{(^ + V X Si) • V„(V«£) + + v x B s n ) ■ V v {^ gi )} 
v M i=1 

"7= E e <+ ^(^ + «x AJ-v^Vwa); 

V " i+j>)i+l 
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£ { 9n ~ 9n 

+ {-d t g n -i ~A=J2( Ei + vx B ^ ■ V *(vWj) + E T (9i,m 

» " i+j=n i+j=n+l 
n-1 

(3.7) + e- 1 ^, /=) + E e'-Mr^, /i) + 1%, £)} 

i=l 

f" 

+ £ e^-^T( gi , f 3 ) - — {E e n + v x Bl) ■ V (VM/n) 

i+j>«+2 v " 

^ n— 1 

- -f= E £ ^ + w x Bi) ■ V v {yfjir n ) + {El +vxB s n ). V„(V7^)} 
--^ E £ i+i " n (^+"x5 i )'V c (VM/i). 

V " i+j>n+l 

Take the inner product of J23J) with {I - Pi} {^^j ■ By (OU|) . the LHS of 
(13.61) is bounded from below by 



5||(I-P0 



fn fn 



e 

On the other hand, for the inner products in the RHS of l|3.6|) . by (|1.4()|l . we have 
{-edtf* - v ■ V x f n ~ dtU-x, ±{I - P x }{fl - /„}) 

< Md t f e n \\ + HV^IU + ||9 t /„_i||} x 1||{I - P 1 }{r n - f n }\\ v ; 

(— ^ E + V X B • ~ £ {l - Pl}{/» - fn}) 



< 



{ + ii^iDdiaii, + iiv^-ii)} x i||{i - f n y\\ u . 



i-\-j—n 
i,3>l 



For the remaining terms, we use Lemma 3.3 in [§] to get the upper bound as 

n— 1 n— 1 

u(E inu + ii^iw • (Edi^n + ii^id + ira + \\ B n\\)} 

i=l i=l 

n— 1 1 

+ (E ll/<H" + X — 1 1^ 1 - P l}{/« - Mil 



i=l 



We therefore conclude that 1 1 1 {I — Pi }{/^— f n } 1 1 w is uniformly bounded. Similarly, 
one can show that i||{I - P 2 }{.9 r e i — 9n}|U is also uniformly bounded. Hence there 
exist fn+i, 9n+i such that 

{I _ Pl }(_____.) ^ fn+u { i_p 2 }(_^____\ ^ gn+1 weakly in ||- 1 1„. 
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Subtracting off L/ n +i, Cg n +i on both sides of (|3.6|) and (|3.7|) . we again isolate the 
zeroth order term. Letting e — ► again, one can deduce (|1.14(l for m = n — 1. 

We now turn to the derivation of various hydrodynamic equations based on (|1.14l) 
and (|1.15() . For the pure Boltzmann case, it is well known that the case m = 1 
yields the celebrated nonlinear incompressible Navier-Stokes equations and m > 1, 
the Navier-Stokes-Fourier system of which derivation can be found in |H| (p. 21-25). 
We take the same path to derive new dissipative PDE's from the rescaled Vlasov- 
Maxwcll-Boltzmann system. We will use numerous results from [§]. First consider 
(|1.14|l when m = — 1. It is equivalent to l|1.21f> and hence we get 

M 2 — 3 _ _ 

(3.8) fi = {pi+v-u 1 + ( — )0i}VM and 9i = o-iVi"- 

When m = 0, electromagnetic field equations become V x B\ = 0, V • B\ = 0, 
V x i?! = 0, and V • i?i = <ri. Thus we may assume B\ = and i?i = V</>i for some 
0i which satisfies A</>i = a\ under the restriction J B\dx — and J E\dx = 0; 
indeed, a nonzero constant B\ does not produce any new macroscopic terms in 
(|1.14ll when m = 1, since 

v x Bi ■ V„( v / /I(jii) = v x Bi ■ V v (aifi) = v x B\ ■ (—a(i)v = 0, 
(_L W x B x ■ V^(Vm/i)i V/ 7 ) = (~i=v x B x ■ \7 v (v ■ uifj,), y/JT) = 0, 

where the last equality is due to the integration by parts in v. Recalling the kernel 
of L and £, by collision invariant property, 

(3.9) {v V x /i,[l,u,H 2 /2]VM) =0 and {v ■ V x gx - E x ■ v^jl, Vm) = 0. 

Plugging l|3.8l) into l|3.9|l . the very first equation gives rise to the incompressibility 
(|1.23|l and Boussinesq relation (|1.24|) . Microscopic part of f 2 and g 2 can be written 
as following by solving l|1.14fl when m — 1 for Lf 2 and Cg 2 : 

{i - Pi}/ 2 = l-'i-v ■ v,/i + r(/i,/i)}, 

(3.10) {I - P 2 }g 2 = C-'i-v ■ V x9l + r(.gi,/i) + ■ 

= ■ (-Vxtri + #i) + oi/i, 

because of the following identity obtained from l|1.13fl and (|1.21|) : 

(3.11) C- 1 T(g 1 J 1 ) = a 1 C- 1 T(^JIJ 1 )=a 1 C- 1 [-Lf 1 +Cf 1 }=a 1 f 1 . 

Notice that {I— Pi}/2 and {I — P2J-.92 are completely determined by already 
known functions. By collision invariant property, we get for m = 1 

(dth + v ■ VJ 2 + -±-E x • V„(Vwi), [M, |^| 2 /2] VM) = 0, 

(3.12) y 

(9t3i + w • V x g 2 + —Ei • V^( v / m/i) - £ 2 • v^M, = 0. 
v M 

The process is similar to the derivation from pure Boltzmann to Navier-Stokes 
except that now we have new terms to deal with; for instance, see P] E] for the 
computation of (d t fi + v ■ \7 x f 2 ,[l,v, \ v\ 2 /2]^/Jl) . Here we pay attention to those 
extras such as electric field related terms. Note that magnetic fields are not involved 
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at this stage. The first equation of 13.12(1 is equivalent to 

dtPi+V ■u 2 + (-^-E 1 • V„(^i),Vm> =0, 
V M 

<9 f ui + ui ■ Vui + Vpi - rjAux + {—Ex ■ V„( v / /Ig i i) ) = 0, 

v M 

^a t {pi + 0i } + | Ul • v0i - -kA^ + • u 2 + (-L^ • v„(Vm0i), ^-Vm) = °- 

Electric field related terms can be taken care of by integrating by parts in v: 
(A= E i ■ V„(v^fli), VM> = ' y V (o-i/i)d» - 0, 

(~^= El ' ^(VWiX'Vm) = Ei- J V v (ain)vdv = -Ei J aifidv = -Exdx, 

(^=Ei ■ V„(VM9i), ^Vm> = E i ■ J Vv(<rili)^-dv = -E x ■ {ax J v^JIdv) = 0. 

Thus the first equation in 1(3.12(1 leads to 1(1.31(1 for m — 2 as well as 1(1.25(1 and 
(|1.28fl . Before going any further, let us integrate 1(3.12(1 over T 3 and then, based on 
the above computation, we get 



— f ui(t,x)dx=— [ 9i(t, x)dx = — [ 
dt J T 3 dt J T 3 dt Jj3 



<Ji(t, x)dx = 



which assures the validity of 1(2.1(1 for m = 1. For the LHS of the second equation 
in 1(3.12(1 . since 

{—Ex ■ V„(Vm/i) ~ E 2 ■ vy/ii, VM> - and (v ■ X7 x P 2 g 2 , ^Jl) = 0, 
V M 

first it is reduced to 

(d t gx + v ■ X7 x g 2 + -^=Ei • V^( v / m/i) - E 2 ■ v^/fl, ^/JT) 

= d t oi + (v ■ V X {I - P 2 }g 2 , y/jM). 

By (13.10(1 . the second term can be written as, 

(v ■ X7 X {I - P 2 }g 2 , ^jl) = {v \7 x {C-\v^I) ■ {-V x ax + E x ) + ax fx}, ^JJ) 
(3.13) = — aAax + aax + V • (<JxUx) 

= -aAui + aux + ux • V<7i, 

where 

a = I £~ 1 (vyfjl) • Vy/Ji dv > 0. 
Hence we obtain 1(1.26(1 and it completes the case m = 1. 
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Next we move onto higher order systems; consider 1(1.141) for m > 2. We shall 
use an induction on to. First by collision invariant property, we get the following: 

(3.14) 

(dtfm + v ■ VcJm+i + -j= ^2 (Ei+vxBi)- V v {y/J2gj), [1, v, 4j-]VK> = 0, 

V " i+j =m + l 

V ^ i+j=m+l 

i,j>l 

= 0, 

where (dtf m + v ■ V x / TO +i, [1, i>, |w| 2 /2]y / /I) was computed for pure Boltzmann case 
in [5] (p. 21-25). Here we will compute electromagnetic field related terms and then 
combine them with the result in HJ. 
Notice that for each i,j > 1, 

H=(^ + v x Bi) ■ V„(Va*&), VM> = 0, 

since each integrand is a perfect derivative in v. Hence the very first equation of 
((3.14(1 yields the (to + l)-th order incompressibility condition: 

(3.15) d t pm + V • Um+l = 0. 

To derive the velocity equation, we look at the next electromagnetic field term: 

i+j=m+l V " 

id>i 

= — J (Ei + v x Bi) gj^/JIdv (by the integration by parts) 

i-\-j=m-\-l 
i,i>l 

= - X! {EiV'E s -(dtEj-x-VxBj)xBi}(fr1H$) 

i-\-j=m-\-l 

i,i>i 

So the second part of the first equation in 1(3.14(1 with (4.8) in [5] is equivalent to 
(3.16) 

dtu m + (v ■ \7 x L^ 1 {-d t f m -i - v ■ VJ m — ^ {Ei + v X. Bi) -V v (yfjigj) 

+ E r (/*« - E • ^ - (^^-i - v x ^) x ^> 

i+j— m+1 i+J— m+1 

i,i>i 

= —V x {Pm+l + #m+l), 

where we have applied (|1.29|) to solve for the microscopic part (I — Pi)/ m +i- Any 
term of which index is lower than to in the above makes a contribution to the 
remainder R^. In particular, all other electromagnetic terms are to be included in 
the remainder except for E m V ■ E\ and E\V ■ E m . For the estimate of the rest, 
we refer [§] (p. 21-25). As splitting u m = PoU m + (I — P )u m , one can readily keep 
track of each term in p. 33(1 , 1(3. 2J1 and (|3.5(l . In addition, to obtain the m-th order 
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Boussinesq relation 1(1. 32J) . we use the pressure p m -i hi EEHt to solve for p m + m . 
Notice that by taking divergence of the (m — l)-th order Navier-Stokes equation 
(|1.33[1 . we have another expression for p m ~i as 

p m -l = A _1 V ■ {-Ui ■ V(P W m -l) - • Vlti + £l<7 m + E m (Ti + ft^-l}- 

For the temperature equation, we need to compute the following electromagnetic 
term: 

i+j=m+l 

= — E] / (Ei + v X Bi) ■ vgj^fjldv (by the integration by parts) 

i+j=m+l J 
i,j>l 

E Ei- J vgj^fjidv (since V X W ■ V = for V, W vectors in R 3 ) 



i-\-j—m-\-l 



5] ^ ■ - v x (by ma) 



m+l 



Then the last part of the first equation in (|3.14|1 by using (4.8) in [H] becomes 
(3.17) 

-d t 6 m + (v ■ V x L~ 1 {-d t f m -i ~ v ■ VJ ra — E (Ei + v X -Vvi^gj) 



i-\-j—m 
i,j>l 

+ E r (/^/i)}'^v / M)+ E ^-(Pm-VxB,) 

rn+1 i-\-j—ni-\-l 

= 9t{p m + m }, 

where (|3.15|) has been used. One can easily deduce p. 35(1 as well as B3.4JI . 
Now we come to the second equation in l(3.14|l . It is equivalent to 

(3.18) d t a m + (v • V a (I - P 2 )<?m+i, Vm> = 0, 

because (-p= X)*+j=m+i(-Ei + v X B 4 ) • V v (y/fifj) - E m+ i ■ Vy/JI,y/fj) = and 

i,j>i 

(v ■ Va;P2.gm+i, y 7 /^) = 0. Solving for the microscopic part (I — 'P^Qm+i by (|1.3U|) 
and plugging it into l|3.18|l . we deduce that 

(3.19) 

d t a m + (v • Vx^H-C 1 - p 2){v ■ V x P 2 g m ) + E m ■ vy/ji + r(0i,Pi/ m ) 

+r(p 23m ,/ 1 )},v^) 

= (« ■ V x £- l {d t (I - P 2 )g m -i + (I - P 3 )(« • V*(I - P 2 )5m) - r(. 9 i, (I - Pi)/ ro ) 
-r((I-P 2 ) 5m ,/ 1 ) + ^ = E (Ei+vxB^-V^fj) 



E r (fi>/j)}> "\/^- 



i+i— m+1 
»,J>2 
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Note that each term in the RHS of 1)3. 19|) essentially has an index of lower than m 
and hence is part of 1)3.3)1 . Terms in the LHS can be computed in the same way as 
in ()3.13)) . Recalling P25m = cr my /JI, we get 

(v ■ V^/T 1 ^ - P2)(« • V x P 2 <7 m )}, ^ffi) = aAa m , 
(v ■ V x £~ 1 {E m ■ v^/JI}, y/JT) = aa m , 

(v • v x £ -1 {r(0i,Pi/m) + r(p 25m , /i)}, y/ji) = v • (a x n m ) + v • (<r m ui), 

where we have used the fact T{ x fJl 1 Pif m ) = CP\f m and T(y/Ji,fi) = £fi which 
follow from the following observation: for any h with Lh = 0, 

Ch = -T(h, y/jj) = r(y/Ji, h), since - Lh = Y{h, y/JJ) + T(yfji, h). 

Therefore, 1)1.34)1 can be derived with the remainder R° n 1)3.3)1 . 

As for field equations, 1)1.15)1 with the definition of V-2,g m immediately leads to 
1)1.36)1 and 1)1.37)1 . To get the compatibility conditions (|1.38l) . take the integral in x 
of field equations in 1)1.15)1 and then the periodic boundary conditions give rise to 

^- [ B m dx = 0, ^- [ E m dx = - { [ vg m+ iy/JI dv dx. 
at j T 3 at j T 3 j T 3 j r3 

Thus J T 3 B m dx — may be assumed for all time. The dynamics of J T3 E m dx is 
more complicated. First note that / J P2<?m^ ^/Jldvdx = 0. By using 1)1.30)1 . we 
have 

d - ( E rn dx = - I [ (I - P 2 )5m+i Vy/JI dvdx 
Jt 3 Jt 3 Jr 3 



dt 



/ / vy/ji C l {-d t g m -i - v ■ V x g m + V" r(g%, f 3 ) 

id>i 

-— ^ (Ei+vxBi)-V v (s/jIfj) + E m -VT/jl}dvdx 

* i-\-j—m 
*,3>1 

E m dx -\- £ m —ij 

J3 



where 



I I d t g m -iC 1 (v x /JI)dvdx 
Jt 3 Jm.3 



( ' ) ' 201 + J2 f f ^(Ei + vxBd'VviyfctfC-^vy/fldvdx. 

Note that l m -\ only contains terms of index lower than m and J E m dx at t = 
for m > 2 can be arbitrarily given. Other terms have vanished after the integration 
because of either a perfect derivative or collision invariant property. 

It remains to verify conservation laws 1)2. 1)1 for m > 2 to finish Lemma. Firstly, 
it is easy to see 

— / p m dx = and — / a m dx = 0, 
at Jj3 at J T 3 

by integrating (|3.15l) and 13.18fl over T 3 . We give the detailed computation for the 
average of the velocity and the similar argument can be applied to the temperature 
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field. Integrate l|3.16J) to get 

(3.21) 4 f u m dx = V / mv • Ej - (dtEj-x - V x £?,) x £,} dx. 

To simplify the RHS, we observe the following: 

/ E t (V ■ Ej) + E 3 {V ■ E{)dx = - / (V x Ei) x Ej + (V x Ej) x £ 4 dx 



dtBi-i x £!j + dtBj-i x _Ej dx, 

T 3 



/ (V x Bj) x B t + (V x Bi) x B 3 dx = (since V • B % = V • B 3 ■. = 0), 
where we have integrated by parts. Thus the RHS of lj3.21|l is equivalent to 

/ d t Bi_i x Ej — d t Ej_i x Bi dx 

i+j=m+1 Jt» 

»,i>i 

and in turn by rearranging indices we deduce 

d f d f 

— / u m dx = — — > Ei x Bj dx. 

dt J .. m dt 4-; 7 T 3 J 

i,i>l 

Moreover, J T3 u m dx — L 3 Pou m dx, since {/ — Po} u m = Vg m for some scalar func- 
tion c/m by the Helmholtz-Hodge decomposition. Therefore, the desired result fol- 
lows. □ 



4. The Diffusive Coefficients 



We now show that we can determine the coefficients ft, fa, ■■■f m ; <?i, 52, ■■■9m and 
Ei, ...E m ; B\, ...B m by giving the initial divergent free part of velocity field u®(x), 
the temperature field 6f (x), the concentration difference field of (x) and the average 
value ei of initial electric field Ef with e\ = 0. 

Proof, (of Theorem \2. 11 ) We use the induction over r. First consider the case r = 1. 
We need to solve the system (|l~23|) - ((l~2^|) with (OSj) and j m = J 0\ = j a\ = 0. 
Once a priori estimate is given, the existence and the uniqueness of a solution 
[ux(t,x),9i(t,x),cri(t,x),(j)i(t,x)] follow from the standard iteration argument via 
the fixed point theorem. Thus it suffices to illustrate a priori energy estimate in 
the subsequence. First we take <9 7 derivatives of H1.25fl - I|1.28fl : 

dtd-yUi + miV9 7 wi + VftyPi — rjAd^ui = — <9 7l uiV9 72 iti + dy(a{V(f>i), 
dtd^ai + uiVd^ai — aA9 7 (Ti -1- ad^ai = — 9 7l itiV9 72 o"i, 
dtd^Oi + uiV9 7 #i — kA<9 7 #i = — <9 7l wiV<9 72 #i, 
Ad 7 (j)i = d~fCTi, V • Mi = 0. 
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The last three summations are over 71+72 =7 with I71I > 1. Let 7 = 0. 
Multiply first three equation in l|4.1J) by [ut, 01, 61], integrate over T 3 and by the 
incompressibility we get 



1 d 
2~dt 



J \ui\ 2 + \<ti\ 2 + \6i\ 2 dx] +77 J \Vu x \ 2 dx + a J |Vcri| 2 da; + K J iVO^dx 
+ |CTi| 2 dx = J criV^i • u\dx. 

Note that the RHS can be absorbed into the LHS assuming M is sufficiently small 
at time t, since 

/ aiVfa ■ Ul dx < (sup| Ul |)||(7i|| • ||V0i|| < CMWaxW 2 , 

J X 

where we have used the Sobolev embedding theorem and the L 2 estimate for the 
poisson equation A<j)i = <j\. Similarly, the energy estimate for I7I < 2 shows 

\j t \J |a y «i| 2 + |a r CT 1 | 2 +|a r fli| 2da ']+'/ / \^d 7 ux\ 2 dx + a J iVd^a^dx 
+k J |V5 7 6»i| 2 dx+a J |9 7 cri| 2 dx 

< supduii + \9i\ + i^Ddi^mii 2 + ||V9 7U1 || 2 + ||a 7 a 1 || 2 + ||va 7 a 1 || 2 + iia^H 2 

X 

+ ||va 7 0i|| 2 + ||v9 7 (/.i|| 2 ). 

Each term in the RHS can be absorbed into the LHS such that 

^(IMIh" + + + |liv^i|| 2 ff2 + f \WiWh + \ llv^H 2 ,, < 0. 

Applying the Poincare inequality, for A = -r min{?7, a, n} > we have 

WuMlm + lki(*)IU> + \\6i(t)\\ H * < Ce- A *{lkill^ + \W°i\\m + 

which immediately implies the existence, uniqueness and exponential decay of 
Ui(t,x), ai(t,x) and 9i(t,x). 

Next we turn to high order derivative cases. We claim that for s > 3, there exists 
a polynomial U s > with ?7 S (0) = such that 

(4.2) |K(t)|| ff . + + |ki(t)|| H . < e- At C/ s (|| U ?||^ +\W\\ H . +W\\h-). 

Separating the case of I'yx] = 1 or I71I = s, and the case of I71I < s — 1, < s — 2, 
we estimate the L 2 norm of the RHS' in (|4.1|l by 

{GdMi^ + iMM + m E (11^x11 + 1^11 + 11^11)+ E IIWill} 

|7|=*+1 |7l=« 

+ C £ {|K|| i? a- 1 +||cT 1 || H a-l+||^|| if a- 1 +||V^|| H a- 1 } 2 , 
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for any small number £ > 0, by an interpolation in the Sobolev Space. We then use 
standard energy estimate for d^ui, d 7 ai and 9 7 #i to get (V • Ui = 0) 

\j t E{ii^ u iii 2 + ii^iii 2 + n^iii 2 } 
i 7 i= s 

+ E {^ll V5 7 u iH 2 + «l|Va 7 0i|| 2 +a||5 7 ai|| 2 +a||Va y ai|| 2 } < 
h\=* 

{c(\\ui\\ H * + iMM+m E ii^iii + HViii + ii^iii + E n v5 7^iii} 2 

h\=s+i h\=s 
+ C c {||«i|| H .-i + IkiH^-! + \\9i\\h-i + HV^H^-i} 4 . 

Note that ||V</>i|| ffa -i < C||o-i|| H .-i. For C(\\ui\\ H 2 + \ \(7i\\ H 2) and £ sufficiently 
small, we deduce from the Poincare inequality and an induction for s — 1 over the 
last lower order term that 

E PWOII + Iiar0i(i)|| + ||a r 1 (t)|| 

l7l<* 

<e- At £ l|07"i(O)ll + 1107*1 (0)|| + ||0 7 0i(O)|| 

|7|<* 

+ f e-^ye- 2 ^{U a - 1 (\\u° 1 \\ H .-i + + \\#l\\H-i)} 2 dT 

Jo 

< e-^U s (\\ Ul (0)\\^ + lki(0)|k- + ll*i(0)||ff.)- 

We thus conclude our claim l|4.2[) . 

Now let us turn to general space-time derivatives d T . Notice that 

Api = V • {(7iV^>i — u\ ■ Viti}. 

We then use repeatedly the equations (|1.25|) , (|1.2(i|) and l|1.28|l to solve for temporal 
derivatives to get 

E {l|0^i(*)H + l|0r*i(*)ll + \\drQx{t)\\} < e- At ^ 2s (|| Ul (0)||^ + lki(0)||^ 
|t|<« 

+ 11^(0)11^.). 

Notice that we need twice many spatial derivatives now for the initial data. Finally, 
since 

_ | w |2 5 _ 
fi = v- ui*/(i+ {— 2^ 1 V / A*j and 9i = ciVfi' 

it follows that for any s > and /3, 

E l|[0^/iW,0^i(«)]||, <CE {H0^iWII + H0r*i(4)|| + ||0^i(4)||} 

|t|<« |t|<s 

< e- A *f/ 2s (|| M 0|| ff2a + + |K||^)- 

Recall that £i = V^>i with A</>i = o\ and Si = 0. Our theorem thus is valid for 
r = 1. 

Assume that fi,fs, --fr, 5ii52,---3r, E\, E2, ...E r , and B\{= 0), B2, ...B r have 
been constructed to satisfy ljl.29|l - l|1.38|l for up to r > 1. We now construct f r +i, 
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g r +i in two steps. 

Step One : Construct the microscopic part {I — Pi}/ r +i and {I— P2}<?r+i 
from the microscopic equation 

In order to solve for {I — Pi}/ r +i and {I — V^gr+i, we need (|3.14ll for rn = 
r — 1 , which are equivalent to fTTK| . J3IIEJ), lpTT7|l and l|3~T3)) for m = r - 1, to be 
true. By the induction hypothesis, the r-th incompressibility condition p. 31(1 and 
the (r — l)-th temperature equation ((1.35|) imply that 1(3.15(1 and 1(3.17(1 hold for 
m = r — 1. And (r — l)-th concentration difference equation 1(1.34(1 leads to ((3.19(1 
for m = i — 1. Finally, applying the r-th order Boussinesq relation 1(1.32(1 into the 
(r — l)-th order Vlasov-Navier-Stokes system ((1.33(1 . we conclude that the pressure 
Pr-i is given by ((3.5(1 for m = i — 1, which also implies that ((3.16(1 is valid for r — 1. 
Therefore, we can solve {I — Pi}/ r +i and {I — P 2 }gr+i from 1(1.14(1 such that 

L{I - Pl}/ r+ l = -dt.fr-1 - V ■ V x / r - ^(£,+»xB,)' ^vW^9j) 



E r (/<>£); 



£{I - P 2 }.g-r+i = -d t g r -i - v ■ \7 x g r 3_ ^ _|_ v x • V v (rfifj) 



V 7 ^ j+j=r 
ij>l 



+E r -v v / JJ+ E r (9t,fj)- 

i-\-j—r-\-l 

We next estimate such microscopic parts. Since two expressions in the above have 
almost same structure, we compute only one term {I — P2}g r +i- The same argu- 
ment can be applied to {I — Pi}/ r +i. Taking d@ derivatives, we get 

(d?£{I - P 2 }g r+1 , d?{I - P 2 } 9 r+i) 

= {-dPQtgr-x - d?{v ■ V x9r }-dP{^- Y,(E* + vx Bi) ■ V„(VM/i)} 

V " i+j =r 



+d?{E r -v^jI}+d? E T (9i,fj), ^{I-P a }Sr+i) 



i+j=r+l 
ij>l 



Applying Lemma 3.3 in [5], we first get {Cr = if [3 = 0) 



(4.3) 



<{-d t d^g r ^ - • V x9r } + d?{E r ■ v^JI}, d?{I - P 2 }g r+ i) 



+ (-^{4= E ^ + w x B ^ ■ v «Vm/,}, - P 2 }g r+ i) 



=(/) + (J/) + (J//). 
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Since |9^{u.v//Z}| is uniformly bounded, by using (|1.4U|) . (J) is bounded by 

{||afi&&.-i|| + \\d?V wgr \\ v + \\dt^ x 9r\\ + \\d T E r \\} x ||6£{I-P a } Sr+1 ||„, 

where \(3±\ = 1. Applying Lemma 3.3 in [S] again, the nonlinear collision term (II) 
is bounded by 

{\\d?9i\\ • \\d?fj\\u + \\d?gi\\u • x \\d?{I - P 2 }g r+1 \\ v . 

To control the electromagnetic terms (III), we introduce the following inequality. 
Claim. There exists C > such that for |r| < s and j3\ < /3, 
+ « x S) • V„(VM/)}, ^5) 
f44) < C E (ll^ll-II^V./ll + II^BH-II^V./ll, 

^ ' ' |?|<max{s,3} 

+ 11^11- II^VID-II^IU. 

The direct result of the claim is that (III) is bounded by 
J2 Eill^'^lll ' IISfV„/ 3 -||„ + II^H • \\dSfiU} x ||#{i - P 2 } 5r+1 |U 

i+j=r \t\ 
i,3>l 

where |r| < max{|r|,3}. From the induction hypothesis i|2.2[l for /i, ...f r ;g±, ,..g r 
and ...E r ;B\, ...B ri consequently, we deduce that for |r| < s, (there are s + 1 
derivatives for f r ,g r ) 

||[^{I-Px}/ r+ i,^{I-P 2 } gr+ i]||„ 

(4 ' 5) < C e - At C/,. + 1 (^{||^|| ff2a + 2+ 4 ( .-„ + ||fi?|| aa . +3+ « (r -0 + ||<7?||^. +3+ 4(r-0». 

»=1 

To complete (|4.5|l . it remains to prove the above claim. 
Proof of Claim: First we note that 

-j=(E + v x B) ■ V*(Vm/) = (E + v x B) ■ V v f - E ■ 
We compute only one term and other terms can be treated similarly. 

(d?{E ■ V„/}, d T g) = J Jd Tl E- (d^ v f)d?gdvdx (where n + r 2 = r) 

<(/ / |9 Tl S-(^V„/)| 2 ^)i||9^|| 



<C( J2 \\drE\\-\\d^f\\) 

|r|<max{s,3} 

At the last step we have applied the Sobolev imbedding theorem. Note that 
min{|ri|, \t2\} + 2 < max{|r|, 3}. Therefore, this finishes the first step. 
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Step Two : Construct the hydrodynamic field Pi/r+i and P2<7r+i i-e. p r +i, 
Ur+it @r+i, Cr+l and the electromagnetic field E r +i and B r +i. 

First of all, we recall that {/ — Po}it r +i = Vg r +i for some scalar function (jy+i 
and hence from (|1.31(l one obtains 

Aq r+ i = -d t p r , 

where {/ — Po}w r +i has average zero. Hence the elliptic L 2 estimates lead that 
J2\ T \<s \ \d T {I - Po}u r +i\ \ is bounded by J2\ T \<s \\ d rdtp r \\, which is bounded by the 
RHS of l|4.5ll by the induction hypothesis. 

We determine p r +i + r +i from the Boussinesq relation l|1.32fl : 

Pr+i + r +l = A" : V • P {- Ul ■ V{P u r } - u r ■ V{P QUl } + R?} 

+ (Mj^,L-l{« • V x f r -l + r(/ X ,/ r )}) - \d X d r -Ur-Ux. 

Notice that by the formula (1.'>.2H and the induction hypothesis, we easily conclude 
HA^V-Podri^H is bounded by the RHS of Hence £ M < S ||S!r{Pr+i+^r+i}|| 
is bounded again by the RHS of Q4.5J1 . 

Finally, to determine the remaining Pou r+ i 1 oy+i, r +i and -Ey+i, B r+ \ we solve 
the (r + l)-th order linear Vlasov-Navier-Stokes-Fourier system (ll.33f) - (ll.38f) . By 
the standard energy method (V • Pou r +i = 0), there is the unique solution [P u r+ i, 
ay+i, r +i> Er+i, Pv+i]- We illustrate the energy estimates in the subsequence. 
First, B r+ i is estimated and other terms will be done all together. From l|1.37|l and 
the Poincare inequality, we get 

\\d T VB r+1 \\ < \\{l-P 2 }d T g r+1 \\ + \\d T d t E r \\, \\B r+1 \\<C\\VB r+1 \\. 

By (|4.5|l and the induction hypothesis, J2\t\<s HcVPv+iH is bounded by the desired 
quantity. 

Next we move onto other terms. Take d T derivatives of l|1.33f) when m = r + 1, 
multiply by d T Pou r +\ and integrate over T 3 to get 

^||3 r P Ur+l|| 2 + 7/||V<9 T P W r+ i|| 2 

= -{d T {ui ■ VP u r +i + P w r +i • Vux}, d T P u r+1 ) 
+ {d T {E\^r+i + ^r+l^l}, d T P u r+1 ) + (d T Rr +1 ,d T P u r+1 ). 

By the Cauchy-Schwartz inequality and Sobolev imbedding theorem, first two terms 
in the RHS can be estimated as follows: 

{d T {ui ■ VP u r+ i + P u r+ i ■ Vwi}, d T P u r +i) 

<l E \\drP u r+1 \\ 2 + c v E H^iii 4 ; 

\t\<s+1 \t\<s+2 

(d T {Eia r+1 + E r+1 ai}, d T P u r+ i) 

<l E \\9rP Q u r+1 \\ 2 + j e ii^Wiii 2 + 1 E ii^+iii 2 

\t\<s+1 |t|<s \t\<s 

+ c ha e ii^iii 4 . 

M<a+1 



VLASOV-MAXWELL-BOLTZMANN DIFFUSIVE LIMIT 27 

In the same fashion, starting with (|1.34|) and (|1.35|l respectively, one gets 

±±\\d T a r+1 \\ 2 + a\\\7d T a r+1 \\ 2 + a\\d T a r+1 \\ 2 

= -{d T {ui ■ Vcr. r+ i + P u r+1 ■ Vcti},<9 t 0y + i) + (d T Rr +1 ,d T (T r+ i), 

<^ E n^Wiii 2 + ^ E n^o^ + iii 2 +c Q Yl H^iH 4 

|t|<s+1 \t\<s \t\<s+2 

+ (d T R^ +1 ,d T a r+ i); 
^\\d T 9 r+1 \\ 2 + K\\Vd T r+1 \\ 2 

= -(d T {ui • V0 r+ i + P u r+1 • V0i}, d T e r+l ) + {d T R e r+1 , d T e r+1 ) 

<~ E n^wiii 2 + ! E n^o^+iii 2 + c K H^iii 4 

|t|<s+1 |t|<s |r|<s+2 

+ (9 T i?^ +1 ,a r r+1 ). 

For the electric field, from l|1.36[) . we obtain 

HflrV^+iH < \\d T d t B r \\ + \\d T a r+1 \\. 

The conservation laws l|2.1[) and (|1.38|) are utilized to handle no derivative terms: 
by Poincare inequality, 



HiW+iH < C\\\7P u r+1 \\ +C\ J P Q u r+1 dx\ 

<c||VPo^ + i|| + c E (ll^ll 2 + M 2 ); 

l<t<r 

lkr+i|| < C\\Va r+1 \\- 

||0 r+ i||<C||V0 r+1 ||+C J2 (IW + M 2 ); 

l<i<r 

\\E r+1 \\<C\\WE r+1 \\ + C\ J E r+1 dx\ 

<C\\WE r+1 \\ + C\e- at J E° +1 dx + e- a{t ~ s h r ds\. 

Based on the above estimates, following the argument in 9 (p. 30-31), by the in- 
duction hypothesis and the Gronwall lemma, one can verify that 112.211 holds for 
m = r + 1 and thus it completes the proof of the theorem. □ 



5. Uniform Spatial Energy Estimate 



In this section, we consider the following model problem: 

d t f £ + - ■ V x f + -U/ e = hi, 
(5.1) e e \ 

d t g £ + - ■ (v x9 £ - y/jiE*) + -£g e = hi, 

e e z 
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coupled with the Maxwell equations 

ed t E £ - V x B £ = - / g £ v^/JIdv + jf, V • B £ = 0, 
(5.2) Jr3 

ed t B e + V x E e = j|, V • E e = / g £ ^JIdv, 

where /if, /i|, jf and jf will be given. We shall establish a uniform space-time 
energy estimate for f £ ,g £ ,E £ and B £ . 

We use a different representation for the hydrodynamic field (fluid) parts Pif £ , 
P 2 g £ in a different way as: 

[Pi/ £ , P29 £ ] = [{a £ (t,x)+b £ (t,x)-v + c £ (t,x)\v\ 2 }^l, d £ (t,x)^JI}. 

Our goal is to estimate a £ (t,x),b £ (t,x),c £ (t,x),d £ (t,x) and E £ (t,x),B £ (t,x) in 
terms of (I - Pi)/ £ and (I - P 2 )g £ . 

We remark that our argument in this section is valid for all e < I. In particular, 
the case of e = 1 yields another proof of the existence of the classical solution for 
the Vlasov-Maxwell-Boltzmann system without using the time derivative. 

Lemma 5.1. Assume f £ 1 g £ ,E £ ,B £ are solutions to the system ]5.1)) . such 
that for all t > 0, 

sen 12/ 



a £ {t,x)dx = 0(e(\\[E £ ,B £ ]\\ 2 (t))) + 0(eA), 

T 3 

b £ (t,x)dx = 0(e(\\[E £ ,B £ }\\ 2 (t))) + 0(eA), 

T 3 

c £ (t,x)dx = 0(e(\\[E £ , B £ ]\\ 2 (t))) + O(eA), 

T 3 



(5.3) 



d £ (t,x)dx = and / B £ (t,x)dx = 0, 

T 3 JT 3 

where we have used the standard big O notation in e and A(t) is a function oft. 
And suppose that \ \E £ \\ 2 + ||B e || 2 < M for some constant M . Then there exists a 
constant C\ > such that 

(5.4) 

\\[d,Pif £ ,d 7 P 2 g £ }\\ 2 <s^ + J \\[d,{l-Pi}r^ 1 {I-P 2 }g £ }\\l 

\y\<N h\<N 

+ C l£ 2 ]T Wd.h^ + e'A 2 + C^iMW 2 + HVjfU 2 + HilH 2 + HVilH 2 ) 

\t\<N-l 

where e < 1 and G{t) is a function of t satisfying 

(5.5) |G(t)|< \\[d,f £ ,d 7 g £ W(t)+ ]T \\[d,E £ ,d 7 B £ ]\\ 2 (t), 

\y\<N h\<N 

and d^hfi is the L 2 projection of [d 7 h\ , 9 7 /i|](t, x, v) on the subspace generated by 

WJi, Viyfjl, ViVjy/Ji, Vi\v\ 2 ^fjl\. 
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Proof. As well illustrated in [Hj , there are two fundamental ingredients in the proof. 
First of all, we use the Local Conservation Laws: By multiplying ^/JL^v^JJi, H 2 a//I 
with the first equation of (|5.1|l , *JJi with the second equation and then integrating 
in u G R 3 , we obtain the following: 

(5.6) 3e 6e 6 2 

6= + I{ Vx a £ + 5V x c £ } = - i<« • V X (I - Pi)/ E , V£) + (ht,v^Ji), 

d e t =-\(v V. T (I - P 2 )<? £ , VM> + v 7 ^)- 

The second ingredient is the study of the Macroscopic Equations: notice that by 
plugging / E = Px/ £ + (I - P^/* and g e = P 2 <f + (I - P 2 ).g £ into O. 

£{af + 6^ ■« + cf|v| 2 } A //I+w- {V x a £ + V K 6 e • u + V^M 2 }^ 

- -{ed t + v ■ V X }(I - Pi)/ e - U(I - Pi)/ £ + eftf , 
e d £ ty /Jl + v {V x d £ - E S }^JI 

= -{ed t + v V a }(I-P 2 )<f 

Fixing i and x, and comparing the coefficients on both sides in front of the k/ju, "fyT^) 



= -{ed t + u • VJ(I - P 2 ). 9 £ - -C(I - P 2 )«f + eh\. 

£ 





\v\ s/JI], we obtain the macroscopic equations as 


(5.7) 


V a c £ = 11 - 


-eh e c , 


(5.8) 


ed t c e + dib\ = If - 


-ehf, 


(5.9) 


dib s + d 3 b\ = lf 3 


+■ ehlj for i ^ j, 


(5.10) 


ed t b\ + d t a e = % 


±eh e M , 


(5.11) 


ed t a e = l e a - 


\-eh e a , 


(5.12) 


ed t d e = l d - 


^eh% 


(5.13) 


V x d £ ~E e = l\- 


-eh%. 



Here by elementary linear algebra, the linear parts l E c , If, l e a , l d and 1% are all 

of the form 

cither (-{edt + v ■ VJ(I - P ± )f , - ~(L(I - P 1 )f e , 
(5.14) I 

or (-{edt + v • VJ(I - P 2 ). 9 £ , C) - - P 2 ). 9 £ , C) 

where £ is a (different) linear combination of the basis [y/ji, v^/JI, ViVjy/Jt, Vi\v\ 2 y/Ji] 
accordingly, while h e c , hf , /if ■, h a , h e d and /if. are defined as (/i e ,C) with same 
choices of £. Notice that 

1107*111 + H^fll + W^KjW + woyhuw + \\d 7 hi\\ + \\d 7 h s d \\ + Wd^KW < c\\d 7 m\. 

The macroscopic equations Ij5.7|l - (|5.11|l have the same structure as the pure 
Boltzmann case. Following the proof of Lemma 6.1 in 9], we can deduce the 
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estimates on \7 x d 7l a e , V x d 11 b e and V x d 7l c e : for |7i| < N — 1, 
^||Va 7l 6 £ || 2 < -j t Jje(I-P 1 )d li r,C ij )-d j d 71 b s dx 

+ ^{||v x 9 71 (i - por ||* + p 71 (i - po/ £ || 2 j 
+ c^Hv.a^u 2 + Hv^^H 2 + lia^p 2 }, 

i||va 7lC £ || 2 < -l^ 3 ( E (i- Pl )a 7i / £ ,Cc) • v,a 7lC £ dx 

+ |{||V^ 71 (I - P0.f H 2 , + ||9 71 (I - Pi)/ e || 2 } 

+ Ce^HVaa^H 2 + \\v x d m c\\ 2 + ||a 71 ^|| 2 }, 

J||Va 7l a £ || 2 <-^{ / ( £ (I- Pl )a 7l / £ ,C) • V x & n a*dx + f ed 7l b E - W^a'dx} 
A at Jj3 Jj3 

+ ^{||V X 9 71 (I - P 1 )r \\l + W&nQ. Pi)f\\l} 
+ Ce 2 {\\d 11 h^ + \\Vd 7l b^}. 

We shall, however, estimate V ' x d 11 d e in the same spirit. 

3 

= d i I a Ti E ] + d n l l + ed n K\ by i5TH) 

= d^d e + ^[<9j<Ve + edjd^hl] (since V • E e = <f ) 

Multiply the above by <9 7l <i e and integrate to get 

(5.15) / |V x <9 7l «f | 2 + |9 7l d £ | 2 dx = V / (<9 7l Ze + £<9 7l /if) • djd 7l d £ dx. 

JT 3 j JT 3 

Recall that 

#7i z e =- (£a t (I-P 2 )a 7l g e ,'yv / / 1 ) 

- (w • V X (I - P 2 )d m g s ,v^JT) - ^<>C(I - P 2 )d 71 <f , v^JT). 

Last two terms are of the desired form. As for the first one, we integrate it by parts 
in the t variable: 

J T3(e<9t(I — P 2 )<9 7l .g e , Vy/JT) ■ djd- yi d £ dx 

= 4 / ( £ ( J ~ p 2)d~ /1 g £ ,v y /jl) ■ djd^dx 
at j T 3 

+ ((I - P 2 )d J d 11 g £ , Vy/JT) ■ ed t d 11 d e dx. 
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Replace edtd 7l d £ with d 7l l^ + e9 7l /i| in the latter integral by the macroscopic 
equation ^5.12[l . Thus 1)5.15)1 leads to the following: 

1 d f 

-(||W|| 2 + ||d 7l( f || 2 ) < -- J^(s(I-P 2 )d 7l g^v^jr ) -djd^dx 
(5.i6) + Sfliv^i - P 2 )(f || 2 + ||d 7l (I - P 2 ). 9 £ || 2 } 



The field estimate, the new ingredient of the proof, comes into play for the need 
of the estimates on no derivative terms a £ , 6 £ and c £ to complete the proof of 
Lemma: from (15.311 we obtain by the Poincare inequality, 

Kll 2 + WW 2 + IKII 2 < l|Va e || 2 + ||W|| 2 + ||Vc £ || 2 
( ' ' +e 2 (\\E £ \\ 2 + \\B £ \\ 2 ) 2 +e 2 A 2 . 

Also it will play an important role not only to close the energy estimates in later 
sections but also to achieve decay rates 1)2.21) and 1)2.8(1 . The intriguing computa- 
tion is carried out in two steps: firstly, the electric field will be estimated via the 
macroscopic equation (|5.13|) and then the Maxwell equations (|5.2|) will give rise to 
the estimate on the magnetic field. 

We start with lfB~T5)) : E e = V ' x d e -1%- eh%. 

Iia^ll 2 ^ / {V x d lt d e -d-tJl-ed^hD-d^dx 
Jt 3 

(5.18) </ (V x d 7l d E -ed 7l h e e )-d m E e dx- d 7l l E e ■ d 7l E £ dx 

Jt 3 ' Jt 3 

The last term needs a special care. 

- / djJl ■ d 11 E e dx = / (e<9 t 9 7l (I - P 2 )g £ , v^/JT) ■ d 11 E £ dx 
Jt 3 Jt 3 

+ f < V -V x d 71 (I-P 2 ) 3 £ + i£(I-P 2 )d 71 ^VM> -d^E'dx 
Jt 3 £ 

= (I) + (II) 

It is easy to see that the second term (JT) is bounded by 
E £ \\ 2 C 



" 71 4 " + ^{HVx(i - P 2 )a 71 /|| 2 + ||(i - p 2 )a 71 /|| 2 }. 

As for the first term (I), we first integrate it by parts in t: 

^- ( (e(I-P 2 )d 7l g e ,v^/Ji) ■ d lt E e dx - [ ((I — P 2 )<9 7l g e , Vy/Jt) ■ ed t d 7l E e dx. 
at J T 3 Jj3 



:S2 
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By using the Maxwell equation H5.2J) to eliminate edtd 11 E e in the second term, we 
get the following: 

-/ ((I - P 2 )d 11 g £ jVy/JT) ■ ed t d 11 E e dx 

JT 3 

= f \((I-P 2 )d m g e ,v^Jl)\ 2 dx + f ((I-P 2 )d J1 g 6 ,v^Jl)-d 7l f 1 dx 
Jt 3 Jt 3 

- [ ((I - P 2 )d ll9 E ,v^Ti) ■ V x Q^BFdx 

< |||(I - P 2 )d 11 f\\l + Iwd^jfW 2 + Q||V,(I - P 2 )d 7l g s \\l +e||9 7lJ B £ || 2 , 

where £ is a fixed small number and we have integrated by parts in x to get the 
last inequality. Hence (|5.18|l leads to 

i||a 7l ^|| 2 + ||a 7l #|| 2 <|^<e(I-P 2 )a 7l3 £ ,^>-3 7l i? £ dx 

(5-i9) +Q{||v :c (i-p 2 )a 71 ^|| 2 + ||(i-p 2 )a 71 ^|| 2 } 

+£ 2 ||5 7l ^|| 2 + i||5 7lJ T|| 2 + e||9 7l i? e || 2 . 

As for the estimate on d 7l B 6 , recall that V x B e = ed t E e + ((I — P 2 )g e ,Vy/]l) —jf. 
Let I72I < N- 2. 

||V x d 12 B £ \\ 2 = / {ed t d- l2 E £ + ({I - P 2 )d 12 g £ ,Vy^T) - d-yjf ) ■ xd 12 B £ dx 
Jt 3 

[ ed l2 E £ ■ V x d l2 B £ dx - [ d l2 E £ ■ V x ed t d j2 B £ dx 

dt Jj3 ~ Jj3 

+ f «(I - P 2 )d 72 g £ ,v^JI) - d l2 j{) ■ V x d l2 B £ dx 
Jt 3 

By l]5.2p . the second term containing edtd l2 B £ in the above can be majorized as 
following: 

/ d l2 E £ ■ (V x V x d l2 E £ - V x d l2 f 2 )dx 
Jt 3 

= f d l2 E £ ■ (V(V • d l2 E £ ) - M l2 E £ ~ V x d l2 f 2 )dx 
Jt 3 

= I d l2 E £ ■ \7d l2 d £ dx + I \\7d^E £ \ 2 dx - / d r ,E £ ■ V x d l2 f 2 dx 
Jt 3 Jt 3 Jt 3 

< f \d l2 E £ \ 2 + \Vd l2 E £ \ 2 dx + \ [ \V8 l2 d"\ 2 + |V x 8 l2 f 2 \ 2 dx 

Jf 3 ^ JT 3 

On the other hand, the last term is bounded by 

Q(||(i-p 2 )a 725 £ || 2 + ||9 72 jTH 2 ) + ei|vxa 72J B £ || 2 

for any small number £. After absorbing ^ 1 1 \7 x d l2 B £ \\ 2 into the LHS, we obtain 
for some 0<C 2 = 1- £<1, 

C* 2 ||V x d l2 B £ \\ 2 <^- [ ed l2 E £ • V x d l2 B £ dx + [ \d l2 E £ \ 2 + \Vd l2 E £ \ 2 dx 

dt Jja Jj3 



(5.20) 



\ f |va 72 ^| 2 + |vx9 72 j|| 2 dx + c c (||(i-p 2 )a 72 /|| 2 + ||a 72 jTII 2 ). 

* Jt 3 
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Letting \ j 2 \ = \ \ - 1 > 0, combine Ij5.19jl with H5.20j) to get for some < C3 < 1, 

C 3 ||<9 7l i; e || 2 <— / (e(I-P 2 )d 11 g £ ,v y /JI) ■ d 7l E £ + £C 2 ed 72 E £ • V x d l2 B £ dx 
at J T 3 

(5.21) +Cai|V x (I - P 2 )d ll9 £ \\l + ||(I - P2)d- ll g £ \\l + ||(I - P2)d 12 g £ \\l} 

+e 2 \\d^hl\\ 2 + \\d 7l jf\\ 2 +^Vxd j2 f 2 \\ 2 . 

Now let us go back to (|5.17|l . After applying l|5.2L)[l and (|5.21|) for I71I = 0, 1 and 
I72I =0 with fixed small £, we have 

C 4 (||£ e || 2 + ||S £ || 2 ) < ^- [ sE £ • V x B £ dx 
at J T 3 

1 n 

(5 22) + j f HI - P2)V, 5 £ , v^JI) ■ VE £ + (e(I - P 2 ) 3 e , vVM) • 

+ £{||V a (I - P 2 )V iC <f|| 2 . + ||(I - P 2 )V^|| 2 + ||(I - P 2 )<f|| 2 } 

+ £ 2 (liv^|| 2 + \\K\\ 2 ) + ||i?|| 2 + ||Vj?|| 2 + HifH 2 + HVjlH 2 , 

where < C4 < 1 and e < 1. This estimate with (|5.17ll finally finishes 15.4|l and 
thus Lemma. □ 



Lemma 5.2. Assume f £ ,g £ ,E £ ,B £ are solutions to the system J^5.1\) . \5.%fy and 

satisfy \5.?J\) . Then there exists constant C\ > 1 such that the following energy 
estimate is valid: 

j t {Ci E iWfof*' a#]|| 2 + |P 7J ^, d.B^j-eSGit)} 
h\<N 

+ s E {^ll[9 7 {I-Pi}/ E ,9 7 {I-P2}g e ]|| 2 + ||[9 7 Pi/ £ ,9 7 P2 5 e ]|| 2 } 

| 7 |<iV 

(5.23) <2d E {{d 1 hi;d 1 f) + {d 1 hl,d 1 g E )}+e 2 5 E H^f|ll 2 
| 7 |<jv M<W-i 

2Ci 
+ — - 

e 



E (lia^Hl^ifll + ll^ll -IMII) 

| 7 |<JV 

£ 2 ^ 2 + ^(llifH 2 + ||Vjf|| 2 + I bill 2 + 1 1 VilH 2 ) 



Proof. The standard <9 7 energy estimates with l|5.1|) and l|5.2[) give rise to 
^{||5 7 / £ || 2 + ||9 7 /|| 2 } + ^{(Ld 7 f £ ,d 7 f £ ) + (Cd,g £ ,d l9 £ )} 

-\{v^ji ■ d 7 E £ ,d ig £ ) = (d^hid.n + (d y h e 2 , d l9 £ ). 

Use (|5.2|) twice to deal with — ^(vy/Jl ■ d 7 E £ , 9 7 g e ): 

-(v^/JJ ■ d 1 E £ ,d 1 g £ ) = - I d- ( E £ ■ {ed t d y E £ - V x d 1 B £ - d 7 jf}da 
£ Jt 3 

r {||5 7 £; e || 2 + ||9 7J B e || 2 }-i / {d^E £ ■ 9 7 if + 8 7 B £ R ■ 8 7 f 2 }da 



e 

ld_ 

2dt 
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But from Lemma fo. II there is a constant C\ > 1 such that 

_5_ 

2e 2 



||[a 7 {i-p 1 }rA{i-p 2 }/]|| 2 , 



| 7 |<JV 

E {8\\faPif E AP29 £ ]\\-eSG(t)}-e 2 6 ]T H^f|ll 2 

1 |t|<JV |-y|<JV— 1 

-e 2 5A 2 ^(Hif || a +||Vif || 2 + IbHI 2 + HVilll 2 ). 
By (|1.20() . multiplying by C\ and collecting terms, we deduce our lemma. □ 



6. The First Order Remainder 

In this section, we prove Theorem 12.21 We study the solution to the kinetic 
equation l|2.3|l . We first establish the spatial energy estimate for d 1 f £ and d 7 g £ . 

Lemma 6.1. Let f £ ,g £ ,E £ ,B £ be solutions to the system \2.ty - \&$ which satisfy 
the following conservation laws: 

[ a E (t,x)dx = | / \E £ {t,x)\ 2 + \B £ {t 1 x)\ 2 dx 1 

Jt 3 2 J T 3 



(6.1) 



b £ (t,x)dx = -e I E £ (t,x) x B £ (t,x)dx, 

T 3 Jf 3 

c e (i, x)dx = — 

T 3 6 



/ \E £ (t,x)\ 2 + \B £ (t,x)\ 2 dx, 

JT 3 



d £ (t,x)dx = and / B £ (t, x)dx = 0. 

T 3 JT 3 

Then for any instant energy functional £jv(/ e , g £ , E £ , B e ) where N > 8, there is a 
constant C > such that 

d 
~dt 



: {C\{ 1 1 [8 7 f £ , d l9 £ ] 1 1 2 + 1 1 [d^, 8 7 B £ ] 1 1 2 } - eSG(t)} 



(6 - 2) + ( 5{^||[a 7 {i-p 1 }r,9 7 {i-p 2 }/]|| 2 + ||[a 7 p 1 r,9 7 p 23 £ ]|| 2 } 

< C{£ N (f s , g s , E £ , B £ )^ + £ N (f £ , g £ , E £ , B £ )}V N {f £ , g £ ). 

Notation: We use £n and T>n instead of £ at(/ £ , g 6 , E £ ,B £ ) and T>N{f £ , g £ ) with- 
out confusion. 

Proof. Note that (f Cj . 1 11 falls into the category of l|5. 3J1 with A — 0. We apply Lemma 
O with 

h\ = -T(f £ ,f £ ) + \e £ ■ vg £ -(E £ + vx B £ ) ■ V v g e , 
e 2 

h £ 2 ee -T(g £ , f £ ) + \e £ ■ vf £ ~(E £ + vx B £ ) ■ V v f £ , 

£ 2 

if = 32 = 0- 
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It suffices to estimate the RHS of (|5.23|l : we only need to estimate e 2 ||9 7 /i^| || 2 and 
{d 1 h\,d 1 f £ ) + (dfh^, d 1 g £ ). Terms related to collision kernel such as ^r(/ e ,/ e ), 
ir(g e ,/ e ) are bounded by C{(£n)% + En}~Dn due to Lemma 7.1 in 9 . Here we 
show that the ones related to the electromagnetic field also have the same bound. 



First we look at the projection part ||<9 7 /if,|| 2 including E £ and B £ : 



(6.3) 



/ [ / d~t{\E e ■ vg £ -(E £ + vx B £ ) ■ V v g £ }(dv] 2 dx 
Jt 3 Jr 3 2 

+ [ [[ d~ ( {)-E £ -vf £ ~(E £ + vx B £ )-X7 v f £ }(dv] 2 dx. 
Jt 3 Jr 3 2 



We will only estimate the first term. Let 7 = 71 + 72 where I7I < N — 1. 

(J) = / [/ {l-d^E 6 -vd 72 g £ -(d^E 6 + vxd 11 B £ )-\7 v d l2 g e }(dv} 2 dx 
Jt 3 Jr 3 2 



It 3 Jr 3 

If |7i| < I72I, we take the sup of d 11 E £ , d 1± B 



{I)< [ \d 7l E £ \ 2 \ [ vd l2 g £ (dv\ 2 dx+ [ \d yi E £ \ 2 \ [ d l2 g £ \7 v (dv\ 2 dx 
Jt 3 Jm. 3 Jt 3 Jr 3 

+ [ \d 11 B £ \ 2 \ [ v x d l2 g £ \7 v (dv\ 2 dx 
Jt 3 Jr 3 

<£| 7l | +2 {/ 1/ vd 72 g £ (dv\ 2 dx+ \ d l2 g £ V v (dv\ 2 dx 
Jt 3 Jr 3 Jt 3 Jr 3 

+ I / v x d l2 g £ V v (dv\ 2 dx} (Sobolev imbedding theorem) 
Jt 3 Jr 3 

< S M+2 { [ \v(\ 2 + \V V (\ 2 + \vV v (\ 2 dv) [ [ \d 72 g £ \ 2 dxdv (Holder ineq) 

Jr 3 Jr 3 Jr 3 



< £| 7l |+2 ^|72| (Note that ( decays exponentially) 
<£ N V N {N>8); 

on the other hand, if I72I < |7i|, we take the sup of d l2 g £ : 



{I) <sup\d l2 g £ \{ [ \d lx E £ \ 2 \ j vQdv\ 2 dx+( \d 7l E £ \ 2 \ [ X7 v (dv\ 2 dx 

+ / |<9 7l £ e | 2 | / !jxV„(*| 2 ^} 
jt 3 Jr 3 



< £>| 72 |+4 £| 7l | 
<£ N V N (N>8). 

Thus we conclude that for I7I < N — 1, 

(6.4) e 2 ||a 7 /if||| 2 < C£ N V N . 



3(5 
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Now wc turn to (d 7 /if , 8 1 f £ ) + (8 7 h £ 2 , 8^g £ ) for | 7 | < N. 

I [ d i{l E£ -vg s - (E £ +v x B £ ) ■ V v g e } ■ 8 1 f £ dvdx 
Jt> Jm 3 2 

+ / / a r {]-E e -vf 8 -{E e +vxB s )-V v f £ }-d ig s dvdx 
Jt 3 Jk 3 2 



(6.5) 



^ I I {-d^E 6 -vd 72 g £ -(d 7l E £ + vxd 11 B £ )-S7 v d l2 g £ }-d 1 f £ dvdx 

+ V / / {-d- n E £ -vd^f £ -{d lt E £ + vxd 11 B £ )-V v d 1 J £ }-d 1 g £ dvdx 

+ / / E £ ■vd 1 f £ 8 1 g £ dvdx~ [ [ {E £ + v x B £ ) ■ \I v {d-yf £ d ig £ )dvdx 
Jt 3 Jm 3 Jt 3 Jr 3 

The worst last (I7I + l)-th derivative term of becomes zero after the integration by 
parts. As we note that v{v) ~ (1 + \v\) and J J vfgdvdx < C||/||i/||fl , ||i/ ) we apply 
the Sobolev imbedding theorem to other terms like the previous argument and then 
we can get the desired bound (£jv) 5 £> n- D 

In order to prove 

Theorem E3 it remains to estimate the velocity derivatives. 

Proof, (of Thex>rem \2.2k ) We notice that for the hydrodynamic part [Pi/ £ , P2<? £ ], 

(6.6) II^PxA^Pa/jll < CWld^ir^-PrfW] 

which has been estimated in Lemma 16.11 It suffices to estimate the remaining 
microscopic part 

[^ 1 {I-P 1 }/ e ,^ 1 {I-P 2 } 5 £ ] 

for | 7 i| + \ fi\ < N, |7i| < N - 1 and /3 > 0. We take 8^ of equations |Q and 
sum over I71I + < N to get 

d t d^(i - p^r + U ■ v x d^(i - p 1 )r + - pi)/ e 
+ {dtd^r + \v ■ v x ^pi/ e + \ ■ Vx^-^n 

= 8^{-T{f £ , f £ ) + l -E £ ■ vg £ -(E £ + vx B £ ) ■ V v g £ }; 
(6.7) 

8 t 8^{l - P 2 )g £ + -v ■ [V^(I - P 2 )g £ - d^{^J)d lx E £ ] + ^<£(I - P 2 )g £ 
+ {8 t 8^P 2 g £ + U ■ V x 8^P 2 g £ + I ^d^v ■ [V x d^g e - 8^(^JI)d 7l E £ ]} 
= d^{-T(g £ J £ ) + \e £ ■ vf -(E £ + vx B £ ) ■ V„/ £ }, 

where \/3i\ = 1. We illustrate the estimate on only {I — P 2 }g £ and the other one 
can be done in the same way. Also we will use many results from UJ. Taking the 
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inner product with d!L {I — P 2 }g e , we get 

\j t {\K{l - P2)g £ || 2 } + - P 2 ) 5 e , <(I - P 2 )g e ) 

~(« • d^(^]I)d 71 E s ,d^(I - P 2 )g s ) + (d t d^P 2 g s + \ ■ \7 x d^P 2 g s 

+ -( P ») dfJl v ■ d^V x d 7l g s ,d^(I-P 2 )g s ) 
(6 8) £ 

-(\^) d ^ v -0^(^)0^, 9£(I-P 2 )<f) 

= i«r(/,/ £ ),^ i (i-p 2 )/) 

+ K{\ EE ■ v -f -(E e +vx B £ ) ■ V„/ e }, fl£(I - P 2 ) 5 £ ). 

Here we only estimate the terms including E £ or B e . We utilize the field estimate 
from the previous section. Other terms can be estimated in the same way as 
presented in [§]• See the proof of Theorem 2.2 there (p. 38-40). The third term on 
the LHS of l|tj.8fl can be taken care of as following: 

- £ (vd^(^JI)d 7l E\d^(l-P 2 )f) 

< C^\\d 11 E e \\ 2 + -^|| 3^(1 - P2)g £ \\l (£ is a small positive fixed number) 

d f 

< C e — / (e(I - P 2 )d 11 g s , Vy/Jj) ■ d lx E e + C 2 ed l2 E £ ■ V x d l2 B e dx 

dt J T 3 

+Q{||V x (I-P 2 )9 7l5 E ||2 + \\(I-P 2 )d^\\l + ||(I- P 2 )d l2 f\\l} +£ N V N 

+ ^||<(I -P 2 )g e \\l (by and El) 

where \-y 2 \ <N-2. One gets the same estimate for (± ( Z)d h v ■ d '^ { y /p)d ll E e , 
fl£ (I- P 2 )g e ), since it is bounded by Q||5 71 ^ e || 2 + - P 2 ). 9 £ || 2 . Define 

G(t) by 

G(t) = Q / ((I - P 2 )d 7l g £ , v^/jj) ■ d^E 6 + C 2 d l2 E £ ■ V x d l2 B £ dx. 

Next, the last term in l|6.8|l will be treated. We decompose f € and g e into Pi/ £ + 
(I - Pi)/ £ and P 2ff £ + (I - P 2 ) 5 e respectively: 

(6.9) 

«{r^ ' vf £ -(E £ + vx B £ ) ■ VJ'}, fl£ (I - P 2 )<f ) 

= ' «Pi/ e - (E £ + v x B £ ) ■ V„P a / £ }, fl* (I - P 2 ).g £ ) 

+ «{ V • u(I - Pi)f e -(E £ + vx B £ ) ■ V„(I - Pi)/ £ }, fl£ (I - P 2 ) 5 e ) 
= W + (^) 

The both terms can be computed in the same spirit as in (|6.5I) . The main concern 
is that the number of derivatives could be N + 1 for the worst case due to the 
Vlasov term [\7 v g e ,\7 v f £ ). As for (I), we recall that the hydrodynamic parts are 
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not affected by the velocity derivative as noted in Ij6.6(l and so we can deduce that it 
is bounded by (£n)^T>n- For the second term (77) , the (7V+ l)-th order derivative 
term together with the similar term stemming from the first equation in (|6.7(l is 
gone by the integration by parts: 

{(E* + vxB°)- V„a£ (I - Pi)g £ , < (I - P 2 )/ £ ) 

+ ((E" +vxB £ ). V o a£ (I - P x )f , fl£ (I - P 2 )g £ ) = 0. 

and hence it is also bounded by (£n)^T>n . 

Combining the above estimates with Theorem 2.2 in |2], we obtain for |7i| < 
TV - 1 and \/3\ + | 7l | = N, 

|{i||[<{I-P 1 }r,<{I-P 2 }/]|| 2 } 

+ ±\\[d^{I-P 1 }r,d^ 1 {I-P 2 }f}\\l 

< c{|p 71 / £ ,s 7l ,f ]|| 2 + ^(||v x (i - P 2 )d 7l <f ii 2 + ||(i - Pijd^wl)} 

+ ejG + C{(£ N )^ +£ N + e 2 }V N . 

Multiplying the above by factor 4 and adding a large multiple K of (|6.2[1 . we get 
the following for < N, \0\ + \ < N and |7i| < N — I, 

^mamw^sflw* + (ii^ii 2 + n9 7J B e n 2 )} - sSG(t)} 

( 6 - 10 ) + 2||[fi£{I - Px}r,<{I- P 2 }.9 £ ]|| 2 - AeG{t)\ + SV N 

< C K {{£ N )? +£ N + e 2 }V N . 

Notice that 

G(t) < C {IP 7 Pi/ £ ,9 7 P 2 5 e ]|| + \\d^\\ + \\d^B e \\} • {II^PxA^P^lH 

| 7 |<JV 

(6.11) +||[a 7 {I - Pi}/ £ ,3 7 {I - P 2 }g s }\\ + \\d 7 E s \\ + \\d^\\}. 
And G(t) has the same bound. We thus can redefine an instant energy by 

£ N (f s ,g e ,E s ,B s ) ee KidiWl&r^dygrW + (P 7 £ e || 2 + ||9 7 B e || 2 )} - e5G(t)} 

(6.12) +2||[<{I - P 1 }f e ) d^{I - P 2 }g e ]\\ 2 - 4eG(i), 

for £ sufficiently small. With such a small e > 0, adjusting the constants in (|6.1U|) 
and applying a standard continuity argument, we thus deduce (|2.6|) by letting £n 
sufficiently small initially. Indeed, our argument is still valid for e = 1, since we 
can choose K, G\ large enough. And that guarantees the global existence of the 
Vlasov-Maxwcll-Boltzmann system. 

Since no estimates for the highest iV-th derivatives of E £ , B e are not available 
as we can see that full (up to A-th) derivatives of g £ are needed in the RHS of 
(I5.21|l to control the (A — l)-th derivative of E e , and hence the dissipation rate is 
weaker than the instant energy, we do not expect exponential decay unlike the pure 
Boltzmann case. To obtain a decay rate (|2.8[) . we use an interpolation argument 
presented in [!£]. We point out that our method is slightly easier because we do 
not have to deal with the time derivative. The key idea is to get a bound for the 
A-th derivatives d 1 E e and d 1 B e with the bound of higher energy. Let k be a 
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positive integer and I7I < N. By an interpolation between Sobolev spaces H N 1 
and H N+k , we get 

\\d,E^ + \\d^\\ 2 < ciii^iiH., + ii^iiS^Hii^ii^ + \m\^ +k }. 

Due to (EH and (|^T|) we have 

ii-^iiS? + ww^ir < c N+k £ N+k (o)H\\E e \\ 2 HN -i + iis £ Hh«-i} 

< C N+k £ N+k (0)i{e^-G + V N }, 
at 

where e sufficiently small and 

(6.13) G(t) = [ {(I-P 2 )d yi g e ,Vy/JJ) ■ d 71 E £ + d^E £ -V x d r2 B £ dx, 

up to constant multiple. In particular, we have |G(i)| < £n- Noting that the other 
part of £jv, i.e. the nonelectromagnetic part is bounded by T>n, a lower bound for 
T>n can be given as following: 

CN^N+kiO)-^^ 1 - ejd < V N . 

It follows that 

(6.14) j t {£ N - eG} + C N , k £ N+k {Qi)-^£^ < 0. 

Letting £jy = £n — eG, we get jj£ n < £n < C£n for some C > 1 and thus (|6.14|) 
becomes 

-£jv + CN^N+kiO)-*?^ 1 < 0. 
After dividing this by £ N h and integrating it over [0, t], we get Q2.8jl . □ 



7. High Order Remainder 

In this section, we establish Theorem 12.31 for n > 2. As noted in the re- 
mainder equations l|1.16f) and l|1.17f) contain singular terms such as zeroth order 
terms or first order terms in e which make it hard to apply the energy estimate 
directly. The difficulty is resolved by introducing new unknowns f R , g £ R , E R and 
B R , which can be obtained from further construction of (n + l)-th coefficients 
f n+ i,g n+ i, E n+ i, B n+ i and (n + 2)-th coefficients f n + 2, ffn+ 2, E n +2 ,B n +2 in the 
diffusive expansion i|1.9fl . The compatibility conditions 11.14|l and (I1.15() will then 
eliminate such severe singularity in e. 

We reformulate the problem. Given /1 ,...,/„ ; gi,...,g n ; E\ ,...,E n ; B\ , . . . , B n de- 
termined by the initial data, we further construct artificial / n +i, g n +\, E n +\, B n+ \ 
and g-n+2, E n+ 2, B n+ 2 by letting m = n + 1 and m = n + 2 in Theorem 12.11 

with constant initial data satisfying 

<n = - J37 E /. E i (*) x B % x ) dx > = °; 
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fs[ E J r3 E?(x)-E<j(x)+B?(x)-B°(x)dx. 



3|T , 

i-\-j—m 



Let us introduce new unknowns /jj, and B|j such that 

Ir= fn~ fn — £fn+l — £ 2 fn+2, E R = E^ — E n — sE n+ \ — £ 2 E n+2 , 

(7.1) 

ff_R = Sri ~ 9n — Efjn+l ~ £ <?n+2, -Br = -B n — £?„ — £-B„+l — £ B n+2 . 



Plugging (|7.1|) into the remainder equations l[1.16jl . (|1.17|) and applying (f 1 . 14|) up 
to m — n as well as ()1.36|) . Ijl.37|l for m = n+ 2 we get for n > 2, 



(7.2) 



(7.3) 



where 



d t r R +l-v x r R +±Lr R =hf, 

dtg R + - £ ■ V X (9 R - V» e r) + ^9% = h l 

ed t E R - V x B R = - [ vg%^JIdv ~ e 3 d t E n+2 , V • B £ R = 0, 
sd t B R + V x E £ R = -s 3 d t B n+2 , V-E%= [ g e R ^JIdv, 

JTR 3 



hi = h 1 (J e R ) + hi(f) + h 1 (E e a ,Btt + 11^,3), 
h% = h 2 (g R ) + h 2 (g) + h 2 (E R ,B R ) + h 2 (E,B), 

and each hi is given as follows: 

(7.4) 

71+2 

M/S) = £"" 2 r(/^, /I) + £ ^ 2 { r (/ fl , A) + m, /I)}, 

i=l 

i-\-j>n-\-3 



s n-l 



h x (E R , B%) = - — (E%+vx B%) ■ V v (g%^JI) 

^ n+2 

-— Y +vxBi)- VMy/jT) + (E R + v x B%) ■ V^^VM)}, 

fci(E,fl)=— ^ ^ £ i+i - n - 1 (^+ U x J B i )- V K ( ft VM), 

* i+j>n+2 

n+2 

ftafcfc) = £- 2 r(^, f R ) + £ e*- 2 ^, /,) + T( gi , /£)}, 

i=l 

h 2 (g) = Y, £ i+j ~ n ~ 2 T(9iJj)-£{dt9n + i+v-V x g n+ 2}-E 2 d t g n+2 , 
i-\-j>n-\-3 

pTl— 1 

h 2 (E% B R ) = —j=-(E% + « x Br) ' V„(/^VM) 

n,+2 

~p E ei_1 {(Bi + vxBi)- V v {f RS fJT) + (E% + v x BJj) • V„(jV/Z)}, 

i= i 
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h 2 (E,B) = --L e^-^iEi + vxB t )- V„(/ jA /m) + eE n+2 ■ Vs /J2. 

* i+j>n+2 

Our goal is to study /J,,, g R , E e R and B e R which are equivalent to f £ , g z n , E^ and 
B e n via (|7.1|l . The procedure closely follows the proofs of the last section. We first 
establish pure spatial energy estimate in Lemma 17. II We write the hydrodynamic 
field parts Pi/jj and ~P 2 g e R as: 

[PliL P 2 g R ] = [{a R (t,x)+b R (t,x)-v + c R (t,x)\v\ 2 }^ d%(t,x)^\. 



Lemma 7.1. Assume that \1.8\) are the solution to the kinetic equation $1.6)) such 
that the following conservation laws hold: 

r ^71 r n+2 „ 

/ a R dx=- \E R \ 2 + \B R \ 2 dx + Y j e l - {E i -E% + B i -B R }dx 
Jt 3 l Jt 3 i=l Jt 3 

+ 1 V s t+ i- n I Ei ■ Ej + Bi ■ Bjdx, 

2 i+j>n+3 
l<ij<n+2 

/ b £ R dx = - e n / E R x Bfjdx - V e i I {E, xB% + E R x B^dx 
Jt 3 Jt 3 " i=1 Jt 3 

(7.5) - E 

r n « n+2 „ 

/ 4dx = - — / \E R \ 2 + \B R \ 2 dx --^e 1 {Ek ■ E% + B t ■ B R }dx 
Jt 3 d Jt 3 6 t~i 

I y e W-n / £. . ^. + Bt ■ B dx, 

b ^kt-L* J-* 3 



6 

i+j>n+3 
l<ij<n+2 



/ d R dx =0 and / B R dx = 0. 

jt 3 ' Jt 3 



Then there is A > so i/ia£ /or any £ > 0, t/iere exists some polynomial U^(0) = 

E {ll^/^^ip + iii^,^]!! 2 }-^)} 

ItI<JV 

(7.6) +(5 E {^ll[5 7 {I-Pi}/^9 7 {I-P2}gail^ + ||[5 7 P 1 /^,9 7 P 23 | s ]|| 2 } 

| 7 |<w 

< e- At C/ c CIIKu^, 6»„., ^}lllir)ie 2 + fjvl 
+ C{(£ N )i +£ N + eU £ (\\\{vi n ,O n ,a n }\\\ N )+zyD t r, 
for e sufficiently small, £ N = £N(f R , g R , E e R , B R )(t) and V N = T> N (f R , g e R ){t). 



Notation: £ N [i] = £ N (f it g it E if Bi)(t). 
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Proof. Note that (|7.5|l falls into the category of l|5.3J) with 

71 + 2 

We apply Lemma HTT1 and Lemma ET2l with jf = —e 3 d t E n+ 2, jf = — £3 dtB n +2 to 
obtain the following, an analogue of (|5.23() : 

(7.7) 

|{Ci X! {ii^/L^ai^+n^i,^]!! 2 }-^)} 

| 7 |<iV 



+ 5 E {^ll[5 7 {I-Pi}/| J ,9 7 {I-P2}.9| ? J||^ + ||[a 7 P 1 /^a 7 P 2 .g| J ]|| 2 } 

| 7 |<JV 

<2C X £ {(8^1,8^%) + (8^1,8^%)} + eH £ ||5 7 ^|| 2 
|t|<JV | 7 |<w-i 

+ 2C 1 e 2 (\\8 J E%\\ • ||W?„ +2 || + n^ii • i^a^+211) 

n+2 

+ e 2 <5 (I I ^] 1 1 2 ) 2 + £ 8 ^(l I [ft^„+2, ftS„+2] 1 1 2 + 1 1 [d t WE n+2 ,d t WB n+2 ] 1 1 2 ) 
=(!) + (77) + (777). 

Now it suffices to estimate (7), (II) and (III) in the above to finish Lemma. Ob- 
serve that by Theorem 12. II 

(7.8) (777) < Ce 2 e- 2Xt U(\\\{u n ,9 n ,a n }\\\ 2 N ). 

As for (77), we use the Cauchy-Schwartz inequality: 

(7Q) (II) <Ce 2 (£ N )ie~ xt U(\\\{u n ,e n ,a n }\\\N) 

<Ce' xt U(\\\{u n ,0 n ,a n }\\\ N ){e 4 + £ N }. 

The estimate for (7) can be done similarly as in the pure Boltzmann case. Electro- 
magnetic related terms do not cause any technical difficulty. Following the proof of 
Lemma 8.1 in ||], we can obtain 

e 2 E l|5 7 ^f||| 2 < e - 2At C/(|||{u„,^ l , ( T„}||| 2 v ){e 2 + £: A r} 
| 7 |<iv-i 

(7.10) +C{e 2 £ N + e 2 e- 2Xt U(\\\{u n ,e n ,a n }\\\ 2 N )}V N ; 

{(d 1 h\ 1 d 1 f R )} + (d 1 hl,d 1 g%) < e- At f/ 5 (|||{u„,^, ( r„}||| 2 v ){e 2 + £:Ar} 

h\<N 

+C{e (£ N ) * + eU ( ( \ \ | {u„ , n , a n } 1 1 1 N ) + S}V N . 
Hence, as adjusting constants and collecting terms, we deduce our lemma. □ 

In order to prove Theorem l2.3l it remains to estimate the velocity derivatives. 
Proof, (of Theorem \2.3[ ) It suffices to estimate just the remaining part 

[^{I-P^/S.a^I-Pahfc]. 
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We follow exactly the same argument in the proof of Theorem 12. 21 Comparing dS t 
derivative of i|7.2[l with (|6.7|l . we notice the only difference is that now and 
dS ± h\ are more complicated. Thus it is enough to estimate only 

«&!,a£{i - Pi}/jj) + (d^hi <{i - p 2 } 9r ). 

Here we illustrate how it works for electromagnetic related terms. The terms con- 
taining hi(f R ), /12(g) or h,2{g e p) can be done in the same way as in [§]. See 
the proof of Theorem 2.3 there (p. 50-51). Recalling l|2.2|) . we first get 

(a£fti(25, S), <{I - Pi}/I) + B), - P 2 }g R ) 

< £ e- At ^(|||{u n ,a„,a4||U)||[^ i {I-P 1 }/-,^ i {I-P 2 } g | s ]|| I , 
<e 2 e~ xt U(\\\{u n ,9 n ,a n }\\\ N )(V N )? 

<e a e- aA *t^(|||{u n ,fl n ,CT n }|||^) + |23 Jv . 

Note that we have also taken into account edl^ 1 (E n+ 2 ■ v^/Jl) in h 2 {E, -B). On the 
other hand, as for h\{E R ,B R ) and h2{E R ,B R ), by splitting f R ,g R into Pi/j^ + 
(I - Pi)/£, P 2 3fl + (I " P2)ff« as in (E2J|, one can get 

{d^M E B, B S R ), <{I - Pi} f R ) + {d?M E R> B s n),d^{I - P 2 }g R ) 
<e(S N )i\\ [< g R , I \* ■ || {I - Pi}f B) <{I - P 2 }.9|j] 1 1, 

n+2 

+ B^2e i - 1 {(£ N [i])i\\ [d^g%, % f R ] I \ v + (€jv)*| I V„<*, < V„/i] 1 1„} 

i=l 

•jll[<{I-Pl}/«»^ 1 {I-P 2 }^]|| I/ 

< e(fiv)3X>Ar + ee- xt U(\\\{u n ,9 n ,a n }\\\ N ){V N + (£ N )*(p N )*} 

< ee-^C/dHiu,,.,^,^}!!!^)^ + {e^)* + e - At C/(| ||{u„, ^ T1 , ^Jl Hat)}^. 

Therefore, as in the proof of Theorem 12.21 with all the above estimates as well 
as the ones in [S] (p. 50-51), we deduce that for any £ > 0, 

j t [K{d{\\ &f R , d l9R ] || 2 + II [d y E R , d 7 B%] || 2 }- sSG(t)} 

(7 n) + 2*|| - Pi}/^,^{I - P 2 }g R \\ 2 ~ teSG{t)] + 5V N 

< e - At C/e(|||{u„,^,a n }||| 2 r ){ £ 2 + £ Ar } 

+ C{ {£ N ) 5 + £ jv + e ^ ( 1 1 1 {u„ , n , ct„ } 1 1 1 N ) + QV N . 
We redefine an equivalent instant energy functional £n as in 

to get 

j f £ N +V N < e' x %(\ ||{u„, n , (7„}| ||^){e 2 + £at} 

+ C{(£ jV ) i +f J v + e^(|||{u n ,0 n ,<7 n }|||jv)+?}X>AT. 

Choose and then fix e and £ so that EJ7g(|||{u n , 6 n , (Jn}|||jv) is sufficiently small. 
We also assume that 

£ N < M 
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sufficiently small such that the coefficient in front of T>n satisfies 
C{(S N )i +£ N + eUz(\\\{u n ,9 n ,a n }\\\ N )+Z} < ~. 

Thus we obtain 

(7.12) j^n + ^-Dn <Ce- xt U(\\\{u n ,6 n ,<j n }\\\ 2 N ){e 2 + £ N }. 
In turn, we have 

^_{ e -CU(\\\{u n ,e n ,<r n }\\\ 2 N ) ft e-^ds^y 

Integrating over t, we deduce 

sup ^(t)<e OT (lll{ u - e ™^>lll^/o~e- Xa ^{^(0) + C' e 2 C7(|||{u„,e n , C T n }|||^)}. 

0<t<oo 

We conclude for e sufficiently small and for some other polynomial U, 

sup fjvW<e a (lll^- CT »>lll^{^(0) + £ 2 ^(|||{u n) e n ,a„}||U)} 

(7.13) Q<t<co 

< M/2. 

A standard continuity argument shows that the hypothesis £/\r(t) < M is valid and 
(|7.13|l is proven. 
Recalling 

/Jl = {/n - /n} + e /n+l + £ 2 /ri+2, = - -^n} + ^E n+1 + E 2 E n+2l 

9r = {ffn ~ .9"} + £ 9n+i + e 2 .9n+2, -Bfj = {B e n - B n } + eB n+ x + e 2 B n+2l 
and by Theorem l2.ll 

£at(£/„ + i + £ 2 ,fn+2, £9n+l + £ 2 9n+2 , ?E n+ i + £ 2 E n+2 ,eB n+ i + £ 2 B n+2 ) 

(7 14) 

<Cs 2 e- 2Xt U(\\\{u n ,6 n ,a n }\\\ N ), 

we thus deduce l|2.12|l for -ft - /„, g E n - g n , E*- E„, B e n - B n . 

To get a decay rate ()2.13(l . we use the interpolation argument as done in the last 
part of the previous section. Let I7I < N and k > 1. Recall that 

iia^ii 2 + H^^H 2 < c{\\e%\\$-i + imw^Himi^ + \\B%w]^ +k }. 

Denote a bound for £n by In, i.e. 

/^ = e ff (lll^»^->lll^{^(0)+e a ^(|||{u n ,fl nj( 7 n }||| Ar )}. 
Due to (523, £23), f7~Tn)l and (ZEfl we have 



<(7 iV+fc (/ J v+fe)t{£^G + X'jv + e- 2A *Z7(|||{u n ,^ ! a n }||| 2 v )7 J v}, 

where e sufficiently small and G(t) is defined in (jCi . 1 Hfl so that \G(t)\ < £n- Noting 
that the other part of En, i.e. the nonelectromagnetic part is bounded by T>n, a 
lower bound for T>n can be given as following: 

CNAlN+kr^S^ ~^ t G-e- 2Xt Um^n,O n ,a n }\\\ 2 N )I N <V N . 
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It follows from (f77T^) that 

(7.15) l{£N-e&}+C N , k (I N+k )-*£y < Ce~ xt I N U(\\\{u n ,8 n ,a n }\\\ N ). 

Letting En = £n — sG, we get ^£n < £n < C£n for some C > 1 and thus (|7.15|) 
becomes 

j^N + CN^lN+k)--*^ < Ce~ M I N t/(|||{u n ,0 n ,<T n }|||^). 

Multiply the above by (1 + v) > we get 

(7.16) |{(1 + l) k £ N (t)} < Q(t) + Ce- Xt (l + ±) k I N U(\\\{u n> e n ,a n }\\\ 2 N ), 
where 

Q(t) = (1 + ^) k - 1 S If {t) - C N , k (I N+k )-i(l + l) k £it(t) 

= (1 + ^) fe - 1 Sv(«){l ~ C-^ fe (W)-i{(l + i)^(t)}i}. 
A; k 

To conclude our theorem, it suffices to verify the following statement: 

Claim. There exists CN+k > such that 

SU P {(1 + y) k £ N (t)} < C N+k I N+k , 

t k 

since the same conclusion is valid for £jy and recalling l|7.1|l . combining with (|T. 14|) . 
we can deduce l|2.13l) . 

Proof of Claim: Let S be the set of t such that Q(t) < 0. Note that Q(t) > 
for sufficiently small M and small t, which implies that S c is nonempty. If S is an 
empty set, namely Q(t) > for all t, we can set Cjq+ k = (CVfc) ^ Let t x £ S. We 
can find t < t x so that t Q £ S c and Q(t) < for to < t < tj.. Integrate 1)7. 16[) from 
<o to <i to get 

(1 + j) k £ N (h) < (1 + j) k £ N (to) + CI N U(\\\{u n ,6 n ,a n }\\\ N ) p e- As (l + ^) k ds 

<(C N , k )- k I N+k + CI N U(\\\{u n ,e n ,a n }\\\ N ) j e~ Xs (l + -) k ds. 

Choosing C N+k = (C Ntk )- k + CU(\\\{u n ,e n ,a n }\\\ N ) f™ e - Xs (l + f) fc ds, we con- 
clude the proof. □ 

Acknowledgments: The author would like to deeply thank Yan Guo for many 
stimulating discussions. 
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